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Executive  Summary 

Validation  and  refinement  of  reaction  mechanisms  and  development  of  new  and  improved 
reduced  mechanisms  describing  the  pyrolysis,  oxidation  and  ignition  of  practical  fuels  requires  a 
reliable  database  of  experimental  combustion  targets.  This  database  must  cover  relevant 
combustion  conditions  and  encompass  realistic  practical  fuels,  fuel  surrogate  mixtures  and 
individual  surrogate  components.  Shock  tubes  and  laser  absorption  diagnostics  are  able  to 
provide  three  important  types  of  experimental  measurements  for  this  database:  ignition  delay 
times,  species  concentration  time-histories  and  elementary  reaction  rate  constants.  These 
experimental  methods  are  the  mainstay  of  this  ARO  research  program  at  Stanford. 

The  primary  scientific  problem  that  this  research  addresses  is  the  need  for  a  very  high-quality 
experimental  database  with  which  to  test  and  refine  detailed  reaction  mechanisms  for  jet  and 
diesel  fuel  surrogates.  Over  the  past  several  years,  we  have  made  several  advances  in  shock 
tube  and  laser  absorption  methodology  that  we  believe  are  strongly  supportive  of  the  desired 
rapid  advancement  of  fuel  research  in  this  area. 

Research  was  performed  in  three  main  areas:  improvements  to  shock  tube  performance, 
ignition  delay  times  and  species  time-histories,  and  high-speed  end-wall  imaging.  A  full  list  of 
all  publications  that  acknowledge  support  from  this  contract  is  attached  to  this  report. 

Summary  overviews  of  each  area  follow  and  have  been  adapted  from  the  key  paper  abstracts. 
These  papers  have  been  included  in  this  final  report. 

Improvements  to  shock  tube  performance 

Several  techniques  have  been  developed  for  obtaining  long,  constant-pressure  test  times  in 
reflected  shock  wave  experiments  in  a  shock  tube,  including  the  use  of  driver  inserts,  driver  gas 
tailoring,  helium  gas  diaphragm  interfaces,  driver  extensions,  and  staged  driver  gas  filling.  These 
techniques  were  explored  under  this  research  program  including  the  most  recent  strategy, 
staged  driver  gas  filling.  Experiments  indicate  that  this  staged  filling  strategy  increases  available 
test  time  by  roughly  20  %  relative  to  single-stage  filling  of  tailored  driver  gas  mixtures,  while 
simultaneously  reducing  the  helium  required  per  shock  by  up  to  85  %.  Using  this  staged  driver 
gas  filling,  in  addition  to  the  other  techniques  listed  above,  post-reflected  shock  test  times  of  up 
to  0.102  s  (102  ms)  at  524  K  and  1.6  atm  have  been  obtained.  Spectroscopically  based 
temperature  measurements  in  non-reactive  mixtures  have  confirmed  that  temperature  and 
pressure  conditions  remain  constant  throughout  the  length  of  these  long  test  duration  trials. 
Finally,  these  strategies  have  been  used  to  measure  low-temperature  n-heptane  ignition  delay 
times.  (Campbell  et  al.  2017). 

Ignition  delay  times  and  species  time-histories 

Ignition  delay  times  were  measured  behind  reflected  shock  waves  in  a  shock  tube  for  a  wide 
variety  of  distillate  and  single  component  fuels  over  a  range  of  temperatures,  pressures  and 
mixtures.  The  fuels  studied  include:  jet  fuels  (JP-5,  JP-8,  and  Jet  A),  rocket  propellants  (RP-2), 
diesel  fuels  (F-76  and  DF-2)  and  gasoline.  A  simple  correlation  was  found  to  describe  the 
ignition  delay  times  for  all  these  fuel/air  experiments  for  equivalence  ratios  near  unity, 
temperatures  from  1000  to  1400  K,  and  pressures  from  6  to  60  atm.  A  simple  correlation  was 
also  found  for  low-fuel-concentration  experiments  diluted  in  argon.  Previously  published 
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ignition  delay  time  data  were  found  to  be  in  good  agreement  with  these  correlations.  For 
several  fuels  studied,  systematic  variations  were  seen  in  the  activation  energy  for  ignition  delay 
time  measurements  with  varying  equivalence  ratio  and  oxygen  concentration.  (Davidson  et  al. 
2017). 

Ignition  delay  times  of  normal  heptane  have  been  measured  at  temperatures  ranging  from  651 
to  823  K  and  at  pressures  between  6.1  and  7.4  atm  at  an  equivalence  ratio  of  0.75  in 
15%02/5%C02/Ar  and  in  15%02/Ar  mixtures  behind  reflected  shock  waves  in  a  shock  tube. 
Time-history  measurements  of  fuel,  OH,  aldehydes  (mostly  CH2O),  CO2,  H2O,  and  temperature 
were  also  measured  under  these  conditions.  These  time-histories  provide  critically  needed 
kinetic  targets  to  test  and  refine  detailed  reaction  mechanisms.  Measurements  were  acquired 
using  a  novel  constrained  reaction  volume  approach,  wherein  a  sliding  gate  valve  confined  the 
reactant  mixture  to  a  region  near  the  end-wall  of  the  shock  tube.  A  staged-driver  gas  filling 
strategy,  combined  with  driver  section  extensions,  driver  inserts,  and  driver  gas  tailoring,  was 
used  to  obtain  constant-pressure  test  times  of  up  to  55  ms,  allowing  observations  of  the 
chemistry  in  the  Negative  Temperature  Coefficient  (NTC)  region.  Experiments  with 
conventional  shock  tube  filling  were  also  performed,  showing  similar  overall  ignition  behavior. 
Comparisons  between  current  data  and  simulations  using  the  Mehl  et  al.  n-Heptane  mechanism 
(2011)  are  provided,  revealing  that  the  mechanism  generally  under-predicts  first-stage  ignition 
delay  times  in  the  NTC  region,  and  that  at  low  temperatures  it  over-predicts  the  extent  of  fuel 
decomposition  during  first  stage  ignition.  (Campbell  et  al.  2015). 

High-speed  endwall  imaging 

A  high-speed  OH*  chemiluminescence  imaging  diagnostic  was  developed  to  image  the 
structure  and  homogeneity  of  combustion  events  behind  reflected  shock  waves  in  the  Stanford 
Constrained  Reaction  Volume  Shock  Tube.  An  intensified  high-repetition-rate  imaging  system 
was  used  to  acquire  images  of  OH*  chemiluminescence  (near  308  nm)  through  a  fused  quartz 
shock  tube  end-wall  window  at  10-33  kHz  during  the  combustion  of  n-heptane  (21  %  02/Ar, 
E.R.=  0.5).  In  general,  the  imaging  technique  enabled  observation  of  the  main  ignition  event  in 
the  core  of  the  shock  tube  that  corresponded  to  typical  markers  of  ignition  (e.g.,  pressure  rise), 
as  well  as  localized  ignition  near  the  wall  that  preceded  the  main  core  ignition  event  for  some 
conditions.  Case  studies  were  performed  to  illustrate  the  utility  of  this  novel  imaging  diagnostic. 
First,  by  comparing  localized  wall  ignition  events  to  the  core  ignition  event,  the  temperature 
homogeneity  of  the  post-reflected  shock  gas  near  the  end-wall  was  estimated  to  be  within  0.5 
%  for  the  test  condition  presented  (T  =  1159K,  P  =  0.25MPa  ).  Second,  the  effect  of  a  recession 
in  the  shock  tube  wall,  created  by  an  observation  window,  on  the  combustion  event  was 
visualized.  Localized  ignition  was  observed  near  the  window,  but  this  disturbance  did  not 
propagate  to  the  core  of  the  shock  tube  before  the  main  ignition  event.  Third,  the  effect  of 
shock  tube  cleanliness  was  investigated  by  conducting  tests  in  which  the  shock  tube  was  not 
cleaned  for  multiple  consecutive  runs.  For  tests  after  no  cleaning  was  performed,  ignition 
events  were  concentrated  in  the  lower  half  of  the  shock  tube.  In  contrast,  when  the  shock  tube 
was  cleaned,  the  ignition  event  was  distributed  around  the  entire  circumference  of  the  shock 
tube;  validating  the  cleaning  procedure.  (Troutman  et  al.  2016). 
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Abstract  Several  techniques  have  been  developed  for 
obtaining  long,  constant-pressure  test  times  in  reflected  shock 
wave  experiments  in  a  shock  tube,  including  the  use  of  driver 
inserts,  driver  gas  tailoring,  helium  gas  diaphragm  interfaces, 
driver  extensions,  and  staged  driver  gas  filling.  These  tech¬ 
niques  are  detailed  here,  including  discussion  on  the  most 
recent  strategy,  staged  driver  gas  filling.  Experiments  indi¬ 
cate  that  this  staged  filling  strategy  increases  available  test 
time  by  roughly  20  %  relative  to  single-stage  filling  of  tai¬ 
lored  driver  gas  mixtures,  while  simultaneously  reducing 
the  helium  required  per  shock  by  up  to  85  %.  This  filling 
scheme  involves  firstly  mixing  a  tailored  helium-nitrogen 
mixture  in  the  driver  section  as  in  conventional  driver  fill¬ 
ing  and,  secondly,  backfilling  a  low-speed-of-sound  gas  such 
as  nitrogen  or  carbon  dioxide  from  a  port  close  to  the 
end  cap  of  the  driver  section.  Using  this  staged  driver  gas 
filling,  in  addition  to  the  other  techniques  listed  above,  post- 
reflected  shock  test  times  of  up  to  0.102  s  (102  ms)  at  524  K 
and  1.6  atm  have  been  obtained.  Spectroscopically  based 
temperature  measurements  in  non-reactive  mixtures  have 
confirmed  that  temperature  and  pressure  conditions  remain 
constant  throughout  the  length  of  these  long  test  duration 
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trials.  Finally,  these  strategies  have  been  used  to  measure 
low-temperature  ^-heptane  ignition  delay  times. 
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1  Introduction 

1.1  Motivation 

The  use  of  shock  tubes  for  studies  of  combustion  chemical 
kinetics  has  enabled  the  precise  determination  of  reaction 
rates,  the  measurement  of  ignition  delay  times,  and  studies 
of  other  combustion  phenomena  for  the  past  60  years.  To  date, 
most  of  these  studies  have  employed  test  times  less  than  about 
15  ms  (e.g.,  [1,2])  and  often  times  less  than  2  ms.  Although 
this  short  duration  has  been  sufficient  to  examine  the  chem¬ 
istry  of  hydrocarbon  fuels  of  practical  significance  at  elevated 
temperatures  (>900  K)  and  high  pressures,  recent  interest  in 
low-temperature  combustion  at  practical  engine  conditions 
(between  600  and  900  K  at  moderate-to-high  pressures  and  a 
range  of  equivalence  ratios),  which  is  characterized  by  much 
longer  timescales,  has  created  a  need  for  greater  test  times  in 
shock  tubes. 

1.2  Previous  research 

Several  researchers  have  contributed  to  the  extension  of 
shock  tube  test  times.  For  instance,  Amadio  et  al.  [3]  and 
Hong  et  al.  [4]  developed  driver  gas  tailoring  relationships  for 
use  with  low-sound-speed  gases  and  convergent  shock  tubes, 
respectively.  Additionally,  Dumitrescu  [5],  Stotz  et  al.  [6], 
and  Hong  et  al.  [7]  modified  the  geometry  of  the  driver  section 
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to  achieve  either  non- attenuating  incident  shock  waves  or  to 
mitigate  pressure  increases  in  the  post-reflected  shock  region. 
Rossmann  et  al.  [8]  developed  a  shock  tunnel  with  a  plunger 
in  the  driver  section,  which  allowed  the  driver  section  length 
to  be  modifled  to  reduce  high-pressure  helium  use.  However, 
this  helium  consumption  reduction  came  at  the  expense  of  a 
decrease  in  test  time  as  well.  Similarly,  Pang  et  al.  [9],  Lam 
et  al.  [10],  and  Gates  et  al.  [11]  used  driver  section  exten¬ 
sions  together  with  driver  gas  tailoring  to  achieve  test  times 
approaching  20,  25,  and  45  ms,  respectively  (see  also  Ref. 
[12]).  Hanson  et  al.  [13],  Zhu  et  al.  [14],  and  Campbell  et  al. 
[15]  demonstrated  a  constrained  reaction  volume  approach 
to  conflne  the  reactive  gas  volume  to  an  area  very  close  to 
the  endwall  of  the  shock  tube,  reducing  the  adverse  effects 
of  reactions  occurring  far  from  the  test  location  at  long  test 
times.  Also,  Campbell  et  al.  [15]  used  a  staged  driver  gas- 
fllling  approach,  in  combination  with  the  other  techniques 
listed  above,  to  achieve  test  times  as  long  as  55  ms. 

Several  groups  have  studied  shock  propagation  phenom¬ 
ena  numerically  or  theoretically.  Hooker  [16]  and  Polachek 
and  Seeger  [17]  examined  shock  interactions  at  gas  inter¬ 
faces;  their  work  led  to  better  understanding  of  driver/driven 
gas  mixing  and  the  need  for  a  lightweight  interface  gas  (i.e., 
helium)  near  the  shock  tube  diaphragm  [4].  Frazier  et  al. 
[18]  numerically  examined  the  effects  of  heat  transfer  on 
post-reflected  shock  average  gas  temperature  and  concluded 
that  heat  transfer  did  not  signiflcantly  reduce  area-averaged 
temperature  for  test  times  up  to  (and  longer  than)  20  ms 
in  8-cm-diameter  or  larger  shock  tubes.  Jacobs  developed  a 
quasi-one-dimensional  Lagrangian  code  for  the  simulation  of 
transient  flow  facilities  such  as  free  piston-driven  shock  tun¬ 
nels  [19].  Also,  a  second-order,  MUSCL-Hancock  method 
Unite  volume  code  with  an  exact  Riemann  solver  has  been 
developed  by  the  University  of  Wisconsin  Shock  Tube  Lab 
(WiSTL)  [20].  This  code  is  able  to  generate  x  —  t  diagrams, 
which  graphically  depict  wave  propagation  in  shock  tubes 
and  help  predict  test  time.  Finally,  Rossmann  et  al.  devel¬ 
oped  an  empirically  informed  computer  routine  capable  of 
generating  x  —  t  diagrams  [8] .  The  routine  used  experimental 
incident  shock  attenuation  data  to  better  model  the  incident 
shock  trajectory  and  also  accounted  for  the  acceleration  of 
the  contact  surface  due  to  boundary  layer  growth. 

1.3  Background 

A  pressure-driven  shock  tube  consists  of  two  parts:  the  driver 
section,  fllled  with  high-pressure  driver  gas  (often  helium), 
separated  by  a  thin  polycarbonate  or  metal  diaphragm  from 
the  driven  section,  which  contains  a  low-pressure  driven 
gas  (often  argon,  although  many  shock  tube  studies,  espe¬ 
cially  those  studying  fuel  ignition,  use  nitrogen  dilution). 
The  experiment  commences  when  the  diaphragm  bursts  due 
to  pressure-induced  strain,  and  an  incident  shock  wave  prop¬ 


agates  in  the  driven  section  toward  the  shock  tube  endwall, 
located  at  the  end  of  the  driven  section.  As  it  moves,  the  shock 
wave  heats  and  compresses  the  gas,  setting  it  in  motion  in 
the  direction  of  the  incident  shock.  When  the  incident  shock 
reaches  the  endwall,  it  reflects  back  toward  the  driver  sec¬ 
tion,  further  compressing  and  heating  the  driven  gas  and  also 
stagnating  the  flow.  The  arrival  of  the  reflected  shock  at  the 
test  location  (typically  1-2  cm  from  the  endwall)  marks  the 
beginning  of  the  test  time.  Meanwhile,  an  expansion  wave 
moves  in  the  driver  section  away  from  the  burst  diaphragm 
toward  the  driver  section’s  end  cap.  After  reaching  the  driver 
section  end  cap,  this  wave  is  reflected  and  moves  toward  the 
driven  section  endwall.  When  the  expansion  wave  arrives  at 
the  test  location,  it  causes  a  decrease  in  pressure,  ending  the 
test  time. 

Anx  —  t  diagram  can  be  used  to  clarify  the  various  inter¬ 
actions  between  the  waves  introduced  above,  as  shown  in 
Fig.  1 .  On  this  plot,  the  use  of  straight  lines  implies  that  ideal 
shock  tube  behavior  is  assumed  (i.e.,  shock  and  expansion 
wave  velocity  attenuation  is  neglected);  this  results  in  a  slight 
underprediction  of  test  time.  The  black  line  represents  the 
incident  and  reflected  shock  waves,  the  blue  and  purple  lines 
denote  the  expansion  head  and  tail  waves,  respectively,  and 
the  red  line  represents  the  contact  surface  between  the  driver 
and  driven  gases.  In  this  flgure,  the  dimensions  shown  repre¬ 
sent  those  of  the  shock  tube  used  in  this  study,  with  its  shortest 
driver  section.  This  shock  tube  has  a  driven  section  length  of 
driven  =  9.73  m  and  a  driver  length  of  L driver  =  3.63  m. 
The  test  time  shown  in  this  flgure  is  denoted  r .  Several  areas 
in  this  flgure  are  denoted  with  numbers,  which  correspond 
to  different  regions  in  space-time  during  the  shock  experi- 


Driver  section  Driven  section 


Fig.  1  Idealized  x  —  t  diagram  for  a  basic  shock  tube  experiment, 
generated  based  on  the  WiSTL  code  [20]  with  thermodynamic  data 
taken  from  Ref.  [21].  Conditions:  Mis  =  1-68  (incident  shock  Mach 
number),  T\  =  =  22  C,  P4  =  1.2  atm,  Pi  =  179  Torr,  P5  = 

2.0  atm,  Ps  =  750  K,  Adriver  =  Adriven  (identical  driver  and  driven 
cross-sectional  areas).  Driver  100  %  He;  driven  100  %  Ar 
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ment.  The  initial,  undisturbed  driven  section  is  termed  Region 
1,  and  that  in  the  driver  section  is  termed  Region  4.  Like¬ 
wise,  the  interval  behind  the  incident  shock  wave  bounded 
first  by  the  contact  surface  and  then  by  the  expansion  wave 
head  is  Region  2.  The  region  behind  the  reflected  shock  is 
Region  5.  Finally,  Region  3  denotes  the  conditions  in  the 
driver  gas  between  the  contact  surface  and  the  expansion  fan, 
and  Region  6  refers  to  the  driver  gas  conditions  following  the 
reflected  expansion  wave. 

When  test  times  longer  than  a  few  milliseconds  are 
attempted,  a  new  phenomenon  occurs  wherein  the  reflected 
shock  interacts  with  the  contact  surface,  sending  either  a 
shock  or  an  expansion  wave  back  to  the  test  location.  This 
causes  either  a  pressure  increase  (shock  wave)  or  a  pressure 
drop  (expansion  wave).  Both  of  these  conditions  are  undesir¬ 
able  and  result  in  the  termination  of  the  test  time  (for  detailed 
V  —  r  diagrams  of  these  cases,  see  Figs.  3  and  4  of  Ref.  [3]).  To 
counteract  this  contact  surface  interaction,  driver  gas  tailor¬ 
ing,  made  possible  by  blending  a  lower- sound- speed  gas  such 
as  nitrogen,  argon,  carbon  dioxide,  or  propane  with  helium, 
is  used  such  that  the  pressures  on  both  sides  of  the  contact 
surface,  following  the  passage  of  the  reflected  shock,  are 
identical.  This  stagnates  the  contact  surface  and  as  a  result  no 
wave  is  reflected  back  to  the  test  section  (see  Refs.  [3,4,22]). 
Ultimately,  the  test  time  is  ended  by  the  arrival  of  the  expan¬ 
sion  wave,  reflected  from  the  driver  section  end  cap,  at  the 
test  location. 

An  idealized  x  —  t  diagram  showing  a  shock  experiment 
at  the  tailored  condition  is  shown  in  Fig.  2.  Region  8  refers  to 
the  time  behind  the  reflected  shock,  bounded  by  the  expan¬ 
sion  wave  head.  Notice  that  the  contact  surface  is  motionless 
(vertical  red  line  between  Regions  5  and  8)  after  interact- 


Driver  section  Driven  section 


Fig.  2  Idealized  x  —  t  diagram  for  a  tailored  shock  tube  experiment, 
generated  based  on  the  WiSTL  code  [20]  with  thermodynamic  data 
taken  from  Ref.  [2 1  ]  and  driver  gas  tailoring  generated  using  the  method¬ 
ology  of  Amadio  et  al.  [3].  Conditions:  Mjs  =  1.68,  Ti  =  T4  =  22  C, 
P4  =  2.0  atm.  Pi  =  179  Torr,  P5  =  2.0  atm,  Ts  =  750  K, 
^driver  =  ^driven  -  Driver  46.2  %  He/53.8  %  CsHg;  driven  100  %  Ar 


Driver  section  Driven  section 


Fig.  3  Idealized  x  —  t  diagram  comparing  two  driver  lengths  (3.63  and 
9.57  m),  generated  based  on  the  WiSTL  code  [20]  with  thermodynamic 
data  taken  from  Ref.  [21]  and  driver  gas  tailoring  generated  using  the 
methodology  of  Amadio  et  al.  [3].  The  conditions  are  identical  to  those 
in  Fig.  2 

ing  with  the  reflected  shock  wave,  but  before  encountering 
the  expansion  wave  head.  Also,  the  test  time  in  this  case  is 
almost  ten  times  longer  than  in  the  untailored  example.  In 
this  figure,  a  driver  mixture  of  helium  and  propane  (see  Ref. 
[3])  was  chosen  to  enlarge  the  area  of  Region  8  at  the  given 
temperature,  making  the  plot  easier  to  read  (tailoring  is  also 
possible  using  helium/nitrogen  at  this  temperature,  although 
this  results  in  a  smaller  Region  8  area). 

An  X  —  t  diagram  can  also  be  used  to  predict  test  time 
as  a  function  of  the  length  of  a  shock  tube’s  driver  section. 
For  instance,  consider  Fig.  3,  an  idealized  x  —  t  plot  showing 
two  different  driver  section  lengths  (3.63  and  9.57  m,  rep¬ 
resenting  two  driver  lengths  available  for  the  shock  tube  in 
this  study).  As  evident,  expanding  the  driver  section  length 
increases  test  time. 

Now  consider  a  fictitious  shock  tube  whose  total  length, 
driven  plus  driver,  is  20  m.  Using  an  idealized  x  —  t  diagram 
allows  the  test  time  of  such  a  shock  tube  to  be  easily  esti¬ 
mated  as  a  function  of  the  driver  length.  We  have  developed 
an  in-house  computer  algorithm  for  this  purpose  that  relies  on 
ideal  shock  relations  taken  from  Refs.  [3,22].  Similar  to  the 
code  developed  by  Rossmann  et  al.  [8],  this  routine  employs 
several  simplifications  that  cause  its  test  time  predictions  to 
deviate  slightly  from  the  experiments.  Despite  this  shortcom¬ 
ing,  we  introduce  it  here  because  it  is  a  useful  tool  for  both 
understanding  trends  observed  in  experiments  and  for  aiding 
shock  tube  design.  Figure  4  shows,  at  one  particular  test  pres¬ 
sure  and  three  test  temperatures,  the  effect  of  increasing  the 
driver  length,  while  decreasing  the  driven  length  to  keep  the 
total  shock  tube  length  constant.  Notice  that  increasing  the 
driver  length,  rather  than  driven  length,  results  in  the  longest 
test  time  at  all  Region  5  temperatures  (T^)  shown.  This  figure 
also  provides  values  predicted  by  the  WiSTL  code  [20] .  The 
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Driven  Length  (m) 


Fig.  4  Approximate  shock  tube  test  time  r  at  the  tailored  condition 
[3]  as  a  function  of  driver  section  length,  for  a  total  shock  tube  length 
of  20  m.  Lines  are  computations  by  the  Stanford  code,  and  individual 
points  are  values  obtained  using  the  WiSTL  code  [20],  with  thermody¬ 
namic  data  taken  from  Ref.  [21].  The  positively  sloped  lines  indicate  that 
increasing  driver  section  length  is  more  effective  at  increasing  test  time 
than  increasing  driven  section  length.  As  discussed  in  Amadio  et  al.  [3], 
tailoring  is  not  possible  at  the  shortest  driver  lengths  at  some  low  temper¬ 
atures.  Conditions:  Ti  =  T/^  =  22  C,  P5  =  2.0  atm,  Adriver  =  Adriven- 
Driver  Yid^2  (tailored);  driven  100  %  Ar 


values  of  both  programs  are  in  good  agreement,  although 
the  simple  assumption  in  the  current  Stanford  calculations 
that  the  expansion  wave  velocity  remains  constant  results 
in  a  longer  test  time  prediction.  It  is  also  important  to  note 
that  both  the  WiSTL  and  Stanford  codes  assume  a  driver-to- 
driven  section  area  ratio  of  4^™  _  t  As  a  consequence, 
these  codes  predict  a  driver  gas  helium  mole  fraction  that  is 
higher  than  that  necessary  in  the  Stanford  tube,  which  has  a 
convergent  driver-driven  design  with  =  2.38  [4,23] 

■^driven 

(see  dimensions  in  Sect.  2).  This  results  in  higher  driver  gas 
sound  speeds  and  corresponds  to  an  underprediction  of  test 
time  relative  to  times  observed  in  the  Stanford  facility  (the 
sound  speeds  and  other  properties  of  several  gases  at  20  C 
are  provided  in  Table  1). 

Plots  can  also  be  generated  to  estimate  the  available  test 
time  as  a  function  of  Region  5  temperature  Ts,  for  different 
driver  section  lengths  and  a  fixed  driven  section  length.  This  is 
shown  in  Fig.  5.  Note  again  the  good  agreement  between  the 
WiSTL  code’s  predictions  and  those  of  the  Stanford  code.  As 
temperature  increases,  available  test  time  decreases,  because 
at  the  tailored  condition  more  helium  is  required  in  the  driver 
gas,  thereby  increasing  the  speed  of  the  expansion  wave 
head.  Conveniently,  above  about  900  K,  the  timescale  for 
kinetic  processes  generally  decreases  with  temperature  as 
well.  This  implies  that  a  shock  tube  with  sufficient  driver 
length  operating  at  the  tailored  condition  is  able  to  observe 


kinetic  processes  over  a  very  large  range  of  temperatures. 
Plots  such  as  Fig.  5  are  useful  in  determining  the  driver  length 
necessary  to  obtain  the  required  test  time. 

In  the  High  Temperature  Gasdynamics  Laboratory  at  Stan¬ 
ford  University,  we  have  developed  numerous  techniques  for 
extending  the  test  time  in  shock  tubes,  including  driver  inserts 
[7],  driver  gas  tailoring  [4],  the  use  of  helium  interface  gas 
in  front  of  the  diaphragm  [4],  driver  extensions  [9-11],  and 
staged  driver  gas  filling  [15].  Our  goal  in  this  work  is  to  use 
these  tools  sequentially  to  achieve  and  even  exceed  the  test 
time  predictions  offered  by  the  x  —  ^  diagrams  shown  above. 


2  Experimental  setup 

As  introduced  above,  the  shock  tube  employed  in  this  study 
was  initially  composed  of  a  9.73-m  driven  section  with  a 
round  internal  diameter  of  11.5  cm,  and  a  driver  section 
composed  of  a  3.63-m  driver  section  with  a  round  internal 
diameter  of  17.8  cm.  Using  driver  extensions,  which  had 
internal  diameters  of  15.4  cm,  the  driver  section  length  was 
extendable  to  a  maximum  of  13.42  m.  The  driver  and  driven 
sections  were  separated  by  a  polycarbonate  diaphragm  with 
a  thickness  of  0.25  mm,  and  a  four-blade  cutter  was  used  to 
decrease  shock- to- shock  variation  in  burst  pressure  and  min¬ 
imize  the  occurrence  of  diaphragm  fragments  released  in  the 
tube.  Driver  inserts,  discussed  later  in  this  article,  were  con¬ 
structed  according  to  the  methodology  of  Hong  et  al.  [7]  to 
counteract  the  non-ideal  pressure  rise  in  the  test  section  asso¬ 
ciated  with  boundary  layer  formation  and  shock  attenuation. 

Diagnostics  were  located  2  cm  away  from  the  endwall  of 
the  driven  section.  Sidewall  pressures  were  monitored  using 
a  Kistler  603B 1  transducer  covered  in  room  temperature  vul¬ 
canizing  silicone  (RTV)  mounted  at  this  location,  and  shock 
speeds  were  measured  using  the  time-of-arrival  technique  via 
a  series  of  seven  PCB  113A26  transducers  (with  a  483B08 
amplifier)  and  extrapolated  to  the  endwall  with  an  accuracy 
of  better  than  ±0.2  %.  Data  from  these  instruments  were 
recorded  using  a  Lab  VIEW  data  acquisition  system  with  a 
measurement  rate  of  10  MHz.  Note  that  test  times  measured 
at  the  2-cm  location  (as  in  these  experiments)  and  on  the 
endwall  (as  in  the  simulations  performed  in  this  study)  are 
essentially  identical.  The  temperature  and  pressure  at  the  end- 
wall  behind  the  incident  and  reflected  shocks  were  computed 
using  1-D  shock  relations  (for  instance,  see  Ref.  [24])  to  an 
accuracy  of  better  than  ±0.6  and  ±1.6%,  respectively,  using 
thermodynamic  data  taken  from  the  Burcat  database  [21]. 

Prior  to  each  experiment,  the  driven  section  was  mechani¬ 
cally  pumped  to  a  pressure  of  less  than  0. 1  Torr,  and  the  driver 
section  was  rough  pumped  to  a  pressure  of  less  than  1  Torr. 
Research-grade  argon,  helium,  nitrogen,  propane,  and  carbon 
dioxide  were  obtained  from  Praxair  Inc.  Spectroscopic-grade 
normal  heptane,  used  in  ignition  experiments,  was  obtained 
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Table  1  Properties  of  several 
gases  at  20  C  using  data 
obtained  from  the  Burcat 

Species 

Molecular 
weight  (g/mol) 

K  =  ^(20C) 

Gas  constant 
(J/kg  K) 

Sound  speed 
(20  C)  (m/s) 

Category 

database  [21] 

H2 

2.02 

1.41 

4124.5 

1303.8 

Diatomic 

He 

4.00 

1.67 

2077.3 

1007.4 

Noble  gas 

Ne 

20.18 

1.67 

412.0 

448.7 

Noble  gas 

CH4 

16.04 

1.31 

518.3 

445.6 

Hydrocarbon 

H2O 

18.02 

1.33 

461.5 

424.1 

Triatomic 

N2 

28.01 

1.40 

296.8 

349.0 

Diatomic 

CO 

28.01 

1.40 

296.8 

349.0 

Diatomic 

NO 

30.01 

1.39 

277.1 

335.5 

Diatomic 

O2 

32.00 

1.40 

259.8 

326.0 

Diatomic 

Ar 

39.95 

1.67 

208.1 

318.9 

Noble  gas 

C2H6 

30.07 

1.19 

276.5 

310.7 

Hydrocarbon 

CO2 

44.01 

1.29 

188.9 

267.4 

Triatomic 

C3H8 

44.10 

1.13 

188.6 

249.8 

Hydrocarbon 

Kr 

83.80 

1.67 

99.2 

220.2 

Noble  gas 

SO2 

64.06 

1.27 

129.8 

219.4 

Triatomic 

Xe 

131.29 

1.67 

63.3 

175.9 

Noble  gas 

C2HF5 

120.02 

1.10 

69.3 

149.2 

Refrigerant 

Rn 

222.00 

1.67 

37.5 

135.3 

Noble  gas 

SF6 

146.06 

1.09 

56.9 

135.2 

Insulator 

C2CI6 

236.74 

1.07 

35.1 

104.7 

Refrigerant 

C4F10 

238.03 

1.05 

34.9 

103.5 

Refrigerant 

lUU 

90 
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Fig.  5  Shock  tube  test  time  r  at  the  tailored  condition  [3]  as  a  function 
of  test  temperature,  for  Ldriven  =  9.73  m  and  several  values  of  Ldriver 
(corresponding  to  available  Stanford  shock  tube  driver  section  config¬ 
urations).  Lines  are  computations  by  the  Stanford  code,  and  individual 
points  are  values  obtained  using  the  WiSTL  code  [20],  with  thermo¬ 
dynamic  data  taken  from  Ref.  [21].  As  discussed  in  Amadio  et  al.  [3], 
tailoring  is  not  possible  at  the  shortest  driver  length  at  some  low  temper¬ 
atures.  Conditions:  Ti  =  T/^  =  22  C,  =  2.0  atm,  Adriver  =  Adriven- 
Driver  Yid^2  (tailored);  driven  100  %  Ar 


from  Sigma  Aldrich  and  was  degassed  prior  to  use.  Binary 
(tailored)  driver  section  mixtures  were  prepared  in  the  driver 
section  by  alternately  filling  two  components  (usually  in  three 
cycles)  and  allowing  the  mixture  to  settle  for  at  least  2  min 
prior  to  bursting  the  diaphragm;  this  technique  has  also  been 
used  by  other  researchers  [3].  Due  to  shock-to- shock  varia¬ 
tions  in  filling  pressures,  the  actual  driver  gas  tailoring  varied 
slightly  between  individual  experiments,  as  the  reported  val¬ 
ues  in  the  plots  contained  in  this  article  suggest.  Finally,  to 
initiate  the  experiment,  the  non-helium  component  (usually 
nitrogen)  was  filled  until  the  diaphragm  ruptured.  Because 
filling  ports  in  our  shock  tube  are  located  near  the  driver  sec¬ 
tion  end  cap  in  all  driver  section  configurations,  this  had  the 
side  effect  of  prolonging  test  time  by  pushing  the  helium- 
containing  mixture  toward  the  diaphragm,  much  like  in  the 
staged  driver  filling  technique  described  in  detail  in  Sect.  4. 


3  Application  of  existing  test  time  extension 
techniques 

In  this  section,  we  will  sequentially  demonstrate  the  use  of 
driver  inserts  [5-7],  driver  gas  tailoring  [4],  helium  interface 
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gas  in  front  of  the  diaphragm  [4] ,  and  driver  extensions  [9-1 1  ] 
to  achieve  a  long  test  time  with  uniform  pressure  throughout 
the  shock  experiment.  Details  of  these  individual  techniques 
may  be  found  in  greater  detail  in  the  referenced  papers. 

3.1  Driver  insert 


Pressure  profiles  from  two  baseline  cases  are  shown  in  Fig.  6: 
one  experiment  using  a  driver  insert  and  the  other  without. 
In  both  experiments,  the  driver  gas  was  100  %  nitrogen  and 
the  driven  gas  was  100  %  argon.  Also,  the  driver  and  driven 
lengths  were  3.63  and  9.73  m,  respectively,  and  the  driver- 
to-driven  cross-sectional  area  ratio  was  =  2.38.  At 

■^driven 

time  zero,  both  pressure  records  show  a  step  change  corre¬ 
sponding  to  the  arrival  of  the  incident  and  reflected  shocks  at 
the  2-cm  test  location  (the  resolution  of  the  plot’s  abscissa  is 
insufficient  to  differentiate  the  two  shocks).  Notice  that  in  the 
shock  taken  without  the  driver  insert,  a  subsequent  normal¬ 
ized  pressure  increase  of  approximately  =  1.8  %/ms 
is  observable  (the  non-uniform  increase  occurring  at  about 
9  ms  in  this  trace  will  be  discussed  later).  Pressure  increases 
such  as  these  are  highly  undesirable,  since  they  isentropically 
increase  the  experimental  temperature  according  to  (1)  [25], 
with  representative  values  at  1000  K  as  shown  in  Fig.  7. 


T 


(1) 


In  this  equation,  Tq  and  Pq  are  the  initial  temperature  and 
pressure  in  Region  5,  and  y  =  ^  is  the  temperature- 


Fig.  6  Comparison  of  shocks  with  and  without  a  driver  insert.  The 
pressure  increase  around  9  ms  is  caused  by  a  driver  gas  tailoring  mis¬ 
match.  Notice  that  before  and  after  this  increase,  the  pressure  in  the 
presence  of  the  driver  insert  is  constant  to  within  ±5  %.  Conditions: 
75  =  808  K,  ^5  =  1.5  atm,  4*^  =  2.38,  Ldriver  =  3.63  m, 

^driven 

7-driven  =  9.73  m.  Driver  100  %  N2;  driven  100  %  Ar 


Fig.  7  Percent  increase  in  test  temperature  for  three  test  gases  at  a  nom¬ 
inal  temperature  of  1000  K,  assuming  isentropic  behavior  according  to 
(1).  Gas  properties  computed  using  thermodynamic  data  obtained  from 
Ref.  [21] 

dependent  heat  capacity  ratio  of  the  driven  gas.  These 
pressure  increases  are  caused  by  factors  such  as  shock  atten¬ 
uation  and  boundary  layer  growth  behind  the  incident  and 
reflected  shock  waves  and  can  be  removed  by  use  of  a  driver 
insert  (see  Refs.  [5-7]).  The  driver  insert  serves  to  gradually 
decrease  the  area  of  the  driver  section  (as  a  function  of  dis¬ 
tance)  from  the  diaphragm  toward  the  driver  section  endwall. 
This  area  change  prematurely  reflects  part  of  the  expansion 
wave  back  to  the  test  section  to  counteract  the  non-ideal 
pressure  increase,  resulting  in  constant  post-reflected  shock 
pressure.  Thus,  in  contrast  to  the  record  without  an  insert,  the 
pressure  trace  with  an  insert  shown  in  Fig.  6  exhibits  a  nor¬ 
malized  pressure  change  of  =  0  %/ms  in  the  time  before 
and  after  the  non-uniform  increase  (from  about  6  to  1 1  ms  in 
this  trace).  Finally,  the  pressure  drop  in  both  pressure  traces 
shown  near  19  ms  represents  the  arrival  of  the  expansion  fan, 
reflected  from  the  end  cap  of  the  driver  section,  at  the  2-cm 
measurement  location. 

3.2  Driver  gas  tailoring 

The  non-ideal  pressure  increase  seen  at  a  time  of  about  9 
ms  in  the  no-insert  pressure  trace  of  Fig.  6  is  caused  by  a 
driver/driven  gas  tailoring  mismatch.  This  is  the  result  of 
the  interaction  between  the  contact  surface  and  the  reflected 
shock  wave,  as  discussed  in  Sect.  1.3.  Typically,  a  pressure 
increase  such  as  this  indicates  that  the  driver  gas  requires 
more  helium  at  the  given  Region  5  temperature  (likewise,  a 
pressure  decrease  indicates  a  need  for  a  higher  mole  fraction 
of  the  low- sound- speed  gas,  in  this  case,  nitrogen).  Figure  8 
compares  the  untailored  shock  experiment  (driver  gas  of 
100  %  N2)  with  a  tailored  experiment,  taken  using  a  driver 
gas  of  31  %  He/N2  (the  shock  tube  geometry  is  identical 
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Fig.  8  Result  of  driver  gas  tailoring.  The  31  %  He/N2  driver  gas  cor¬ 
responds  to  the  tailored  condition  at  Ts  ~  800  K  in  a  shock  tube  which 
has  a  driver/driven  area  ratio  =  2.38  [4].  Shock  tube  dimensions 

^driven 

are  identical  to  those  given  in  Fig.  6,  and  a  driver  insert  was  used  to 
mitigate  non-ideal  pressure  increases 

to  that  of  Fig.  6).  At  the  tailored  condition,  the  pressure  in 
Regions  5  and  8  in  the  shock  tube  are  matched  [3,4,22]  (see 
Fig.  2),  such  that  neither  an  expansion  wave  nor  a  shock 
wave  is  transmitted  back  to  the  test  section.  Not  surprisingly, 
a  sustained  pressure  increase  is  not  visible,  indicating  that  the 
reflected  shock  wave  did  not  reflect  from  the  driver/driven  gas 
contact  surface.  However,  a  pressure  hump  (contact  surface 
bump)  can  still  be  seen  near  about  10  ms.  Finally,  the  test 
time  at  the  tailored  condition  is  slightly  reduced,  indicated 
by  the  arrival  of  the  expansion  wave  at  the  2-cm  measure¬ 
ment  location  about  2  ms  earlier.  This  is  the  result  of  helium 
increasing  the  sound  speed  of  the  driver  gas,  allowing  the 
expansion  wave  head  to  move  more  quickly. 

3.3  Helium  interface  gas 

In  an  ideal  shock  tube,  the  contact  surface  between  the  driver 
and  driven  gases  would  have  zero  thickness  and  no  mixing 
would  occur  between  these  fluids.  However,  for  a  variety  of 
reasons  related  to  the  physics  of  diaphragm  rupture,  some 
driver/driven  gas  mixing  does  occur,  which  results  in  a  con¬ 
tact  mixing  zone  of  finite  length.  The  contact  surface  bump 
observable  in  Fig.  8  is  the  result  of  the  reflected  shock  inter¬ 
acting  with  this  mixing  zone,  as  described  in  detail  in  Refs. 
[4,16,17].  Use  of  a  lightweight  gas  such  as  helium  as  a 
buffer  between  the  driver  and  driven  gases  greatly  reduces 
the  influence  of  this  interaction.  In  these  experiments,  helium 
gas  is  introduced  through  a  port  immediately  in  front  of  the 
diaphragm  by  first  filling  the  tube  with  the  driven  gas  mixture 
(argon)  to  a  pressure  very  close  to  the  initial  pressure  (called 
Pi )  and  then  by  filling  helium  until  the  pressure  is  Pi .  Fig¬ 


Fig.  9  Effect  of  helium  interface  gas.  Shock  tube  geometry  is  identical 
to  that  provided  in  Fig.  6,  and  a  driver  insert  was  used  to  mitigate  non¬ 
ideal  pressure  increases 

ure  9  shows  the  impact  of  a  helium  interface  gas  which  is 

2.4  %  of  Pi  (corresponding  to  a  length  in  the  driven  sec¬ 
tion  of  about  23  cm)  on  the  contact  surface  bump.  Helium 
percentages  ranging  from  about  2  to  10  %  give  the  best  per¬ 
formance  (i.e.,  elimination  of  the  contact  surface  bump)  in 
these  experiments.  Note  that  other  researchers  using  differ¬ 
ent  shock  tubes  have  reported  values  as  high  as  20  %  [4] ;  this 
large  range  is  likely  due  to  differences  in  shock  tube  geom¬ 
etry  and  diaphragm  material  (polycarbonate  versus  metal). 
Finally,  observe  that  the  fiat  =  0  %/ms  test  time  in  the 
resulting  shock  is  about  16.5  ms. 

3.4  Driver  extensions 

As  introduced  in  Fig.  3,  lengthening  the  driver  section  can 
greatly  increase  test  time  (see  Refs.  [9-11]).  Figure  10  shows 
the  result  of  increasing  the  driver  section  length  from  3.63  to 
7.74  m,  and  then  to  9.57  m.  True  to  expectations,  increasing 
driver  section  length  delays  the  arrival  of  the  expansion  wave 
at  the  endwall;  test  time  is  plotted  as  a  function  of  length 
in  Fig.  11.  Also  shown  in  Fig.  11  are  test  time  predictions 
made  by  the  Stanford  code.  Notice  that  the  measured  test 
times  are  longer  than  the  predicted  times,  which  is  likely  due 
to  two  factors.  First,  bursting  the  diaphragm  with  nitrogen 
tends  to  increase  test  time  as  discussed  in  Sect.  2;  the  model 
assumed  that  the  entire  driver  section  was  filled  with  the 
tailored  gas  mixture.  Second,  the  shock  tube  employed  has 
a  convergent  design  with  =  2.38,  which  allows  less 

-^driven 

driver  gas  helium  to  be  used  as  compared  to  =  1 

(as  assumed  by  the  models)  to  achieve  the  same  Region  5 
temperature,  thus  decreasing  the  driver  gas  sound  speed  and 
delaying  the  arrival  of  the  expansion  fan  at  the  test  section 
[4,23]. 
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Fig.  10  Comparison  of  pressure  profiles  for  =  2.38,  Ldriven  = 
9.73  m,  with  Ldriver  values  of  3.63,  7.74,  and  9.57  m.  The  driven  gas  is 
100  %  Ar,  and  the  driver  gases  are  tailored  He/N2  mixtures.  A  driver 
insert  was  used  to  mitigate  non-ideal  pressure  increases 


Fig.  11  Test  time  r  as  a  function  of  driver  length,  for  L driven  =  9.73  m 
at  =  800  K  and  =  1.9  atm.  The  driven  gas  is  100  %  Ar,  and  the 
driver  gases  are  tailored  He/N2  mixtures.  These  experiments  employed 
a  driver  insert.  The  line  is  a  computation  by  the  Stanford  code  using 
thermodynamic  data  taken  from  Ref.  [21] 


4  Staged  driver  gas  filling 
4.1  Theory 

Returning  to  Fig.  2,  it  is  evident  that  the  expansion  wave  head 
originates  at  the  diaphragm,  propagates  toward  the  driver  sec¬ 
tion  end  cap,  and  then  is  reflected.  While  undergoing  this 
process,  this  wave  travels  at  the  local  sound  speed  (relative 
to  the  gas)  of  the  tailored  driver  gas  mixture  [22].  Thus,  the 
use  of  helium,  which  has  a  high  sound  speed  (see  Table  1), 
in  the  driver  gas  mixture  impedes  efforts  to  achieve  long 


test  times,  because  it  causes  the  expansion  wave  to  travel 
more  quickly  and  reach  the  test  section  sooner.  Moreover,  the 
price  of  helium  has  skyrocketed  in  recent  years  [26],  mak¬ 
ing  it  costly  to  All  the  entire  driver  section  (especially  once 
this  section  is  elongated  to  expand  the  test  time).  However, 
helium  is  a  desirable  driver  gas  because  of  its  propensity  to 
create  strong  shocks  with  moderate  driver  section  pressures 
(the  only  other  alternative  considered  is  hydrogen  [27-29], 
which  is  not  commonly  used  in  the  USA  because  of  its  explo¬ 
sive  potential).  To  accommodate  the  need  for  strong  shocks 
(obtained  by  having  the  tailored  helium/low- sound- speed  gas 
adjacent  to  the  diaphragm)  while  achieving  the  longest  pos¬ 
sible  test  time  (obtained  by  Ailing  a  majority  of  the  driver 
section  with  a  low-sound-speed  gas),  we  have  developed  a 
staged  driver  gas-fllling  technique. 

Note  that  a  lower-sound-speed  gas  (such  as  nitrogen,  kryp¬ 
ton,  or  xenon)  in  the  driven  section  could  also  be  used  to 
increase  test  time.  However,  nitrogen  has  a  finite  vibrational 
relaxation  time  behind  shock  waves  that  can  affect  data  analy¬ 
sis,  and  krypton  and  xenon  are  prohibitively  expensive;  for 
these  reasons,  argon  is  often  the  preferred  inert  driven  bath 
gas  for  shock  tube  experiments. 

4.2  Technique 

In  a  conventional  shock  tube  experiment  that  employs  driver 
gas  tailoring  such  as  those  introduced  above,  the  driver  sec¬ 
tion  is  filled  with  a  binary  mixture,  often  consisting  of  helium 
and  nitrogen,  until  the  diaphragm  ruptures.  The  staged  Ailing 
technique  involves  one  additional  step  and  is  performed  as 
follows  (see  also  Fig.  12).  First,  the  driver  section  is  Ailed 
to  an  intermediate  pressure  P^a  as  in  the  conventional  AA- 
ing  technique  with  a  tailored  binary  mixture  of  helium  and 
nitrogen.  Second,  the  driver  section  is  Ailed  until  diaphragm 
rupture  (said  to  occur  at  pressure  F4)  using  a  port  near  the 
driver  section  end  cap  (far  from  the  diaphragm)  with  a  third, 
low- sound- speed  (backflll)  gas.  As  this  low- sound- speed  gas 
Alls,  it  pushes  the  binary  mixture  toward  the  diaphragm.  A 
mixing  zone  also  develops  between  the  backflll  and  binary 
gases  during  this  process.  Thus,  the  ratio  ^  governs  the 
extent  to  which  the  driver  section  is  flUed  with  the  binary 
mixture  according  to  (2)  (where  Lbinary  is  the  idealized 
length  of  binary  gas  measured  from  the  diaphragm).  Note 
that  this  equation  ignores  the  presence  of  the  driver  insert, 
which  slightly  decreases  the  driver  section  area  as  discussed 
above.  The  Lbinary  values  provided  in  this  paper  do,  however, 
account  for  the  slight  change  in  diameter  of  our  driver  section 
as  it  changes  from  17.8  cm  over  the  flrst  3.63  m  to  15.4  cm 
over  the  next  9.79  m  (the  total  maximum  driver  length  is 
13.42  m). 


PAu  Lbinary 

Pa  L  driver 


(2) 
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Fig.  12  Staged  driver  filling 
procedure  [in  this  case,  the 
tailored  (binary)  mixture  is 
He/N2].  The  driver  section  is 
first  filled  to  pressure  P/^a  with 
the  binary  mixture  and  then 
filled  from  a  port  near  the  driver 
end  cap  with  the 
low- sound- speed  gas  until 
diaphragm  rupture  at  pressure 
Pa 


^binary 


^4a  ^binary 


P  I 

^  4  '-driver 


^driver 


An  idealized  x  —  t  diagram  comparing  staged  driver  filling 
to  conventional  driver  filling  is  shown  in  Fig.  13.  In  this  plot, 
^binary  —  3.63  m,  Z^driver  —  9.57  m,  and  Z^driven  —  9.73  m. 
For  clarity,  regions  labeled  “a”  correspond  to  the  backfilled 
gas  and  those  labeled  “b”  correspond  to  the  binary  (tai¬ 
lored)  gas.  In  the  staged  driver  gas  shock  (shown  in  solid 
lines),  within  the  driver  section,  notice  that  the  lines  denot¬ 
ing  the  expansion  fan  wave  head  propagating  in  the  backfilled 
(low  sound  speed)  gas  have  a  greater  slope  (time  divided 
by  distance)  than  those  showing  the  conventional  driver  fill¬ 
ing  shock  expansion  fan  wave  head  (shown  in  dashed  lines). 
This  reflects  the  fact  that  in  the  staged  driver  gas  shock,  the 
expansion  fan  travels  more  slowly.  Also  notice  that  within  the 
tailored  gas,  the  initial  slope  of  the  expansion  head  and  tail 
lines  are  low,  reflecting  the  fact  that  the  driver  gas  immedi¬ 
ately  adjacent  to  the  diaphragm  has  a  high  sound  speed.  This 
high  sound  speed  allows  a  strong  shock  wave  to  be  gener¬ 
ated.  The  net  effect  of  this  staged  driver  gas-filling  procedure 
is  that  a  strong  shock  wave  is  produced  via  the  binary  mix¬ 
ture,  but  the  expansion  wave  head  is  delayed  in  reaching  the 
test  location  at  the  end  of  the  driven  section  because  it  must 
first  pass  through  the  low-sound-speed  backfilled  gas  in  the 
driver  section.  Moreover,  helium  gas  is  only  needed  in  a  short 
length  of  the  driver  section  next  to  the  diaphragm;  the  rest  of 
the  driver  section  can  be  filled  with  a  much  less  expensive 
alternative  such  as  nitrogen,  argon,  or  carbon  dioxide  [3]. 

In  the  current  experiments,  we  found  no  evidence  that 
abrupt  expansion  or  shock  waves  are  reflected  back  to 
the  test  section  as  the  reflected  shock  interacts  with  the 
backfill/binary  gas  contact  surface.  In  contrast,  a  tailoring 
mismatch  at  the  tailored  driver/driven  gas  interface  gener¬ 
ates  a  distinct  pressure  change  in  the  test  section,  such  as  that 
demonstrated  in  Fig.  8.  One  possible  explanation  for  this 
discrepancy  is  that  the  mixing  zone  which  develops  between 
the  backfill  and  binary  gases  as  the  backfill  gas  flows  into  the 
driver  section  is  much  more  extensive  than  that  which  devel¬ 
ops  between  the  binary  and  driven  gases  as  the  diaphragm 
bursts.  Because  this  mixing  zone  is  longer,  the  transition  in 
gas  properties  as  seen  by  the  reflected  shock  entering  the 


Driver  section 


Driven  section 


Fig.  13  Idealized  x  —  t  diagram  comparing  staged  driver  filling  (solid 
lines)  to  conventional  driver  filling  (dashes),  with  Zbinary  =  3.63  m, 
Zdriver  =  9.57  m,  and  Zdriven  =  9.73  m.  Notice  that  the  expansion 
wave  travels  more  slowly  in  the  low- sound- speed  backfill  gas,  delaying 
the  arrival  of  the  expansion  wave  at  the  endwall.  Conditions  are  identical 
to  those  in  Fig.  2 


backfill/binary  mixing  zone  is  gradual,  and  a  reflected  expan¬ 
sion  or  shock  wave  is  not  generated  and  sent  back  to  the  test 
section.  For  this  reason.  Fig.  13  is  drawn  with  a  vertical  con¬ 
tact  surface  line  between  Regions  8a  and  8b,  meaning  that 
the  backfill/binary  gas  contact  surface  is  stagnated  by  the 
reflected  shock.  Interestingly,  no  contact  surface  bump,  as  in 
that  shown  in  Fig.  9,  has  been  observed  as  the  result  of  the 
backfill/binary  gas  contact  mixing  zone;  we  attribute  this  to 
the  similarity  in  properties  between  these  gases  (in  the  case 
of  the  binary/driven  gas  contact  mixing  zone,  gas  properties 
can  vary  significantly,  and  the  mixing  of  these  gases  does  not 
necessarily  result  in  linearly  varying  conditions).  Finally,  in 
our  staged  driver  gas-filling  experiments,  we  have  observed 
significantly  lower  non-ideal  pressure  change  when  the 
expansion  wave  propagates  through  the  backfill  gas  as  com¬ 
pared  to  measured  during  conventional  driver  filling  (no 
backfill).  We  attribute  this  to  the  interaction  of  the  expansion 
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Fig.  14  Comparison  of  conventional  to  staged  driver  gas  filling.  The 
use  of  staged  driver  gas  filling  increased  test  time  by  roughly  8  ms  (22  %) 
and  used  85  %  less  helium,  compared  to  its  conventional  counterpart. 
Conditions:  =  2.38,  Ldriver  =  7.74  m,  Ldriven  =  9.73  m.  Driver 

^driven 

30  %  He/N2;  driven  100  %  Ar;  backfill  100  %  N2.  A  driver  insert  was 
used  to  mitigate  non-ideal  pressure  increases 

wave  with  the  mixing  zone  between  the  backfill  and  tailored 
gases;  however,  further  experimentation  is  necessary  to  fully 
understand  this  phenomenon. 

4.3  Staged  vs.  conventional  filling 

Figure  14  offers  a  comparison  between  the  pressure  traces  of 
a  shock  with  conventional  driver  filling  (Lbinary  =  Ldriver  = 
7.74  m)  and  a  shock  with  Lbinary  =  1.1m,  where  the  backfill 
gas  was  pure  nitrogen.  The  conventional  filling  shock  was 
accomplished  by  filling  a  pre-mixed  30  %  He/N2  mixture 
into  the  driver  section.  By  observing  the  delayed  arrival  of 
the  expansion  wave,  it  becomes  evident  that  staged  driver 
gas  filling  expanded  the  available  test  time  by  approximately 
8  ms,  or  about  22  %.  Moreover,  the  staged  driver  filling  shock 
required  about  85  %  less  helium  than  its  counterpart. 

4.4  Variation  of  binary  length 

Figure  15  shows  the  relationship  between  test  time  and 
Lbinary  tailored  He/N2  binary  gas  and  pure  nitrogen  back¬ 
fill  gas.  As  Lbinary  decreases,  a  greater  proportion  of  the  driver 
section  is  filled  with  the  low-sound-speed  backfill  gas  and  test 
time  increases.  Note  that  experimental  Lbinary  values  greater 
than  about  5  m  (except  the  case  of  filling  the  entire  driven 
section  with  a  pre-mixed  binary  mixture)  were  unachievable, 
because  the  corresponding  intermediate  pressure  P4a  was 
too  close  to  the  diaphragm  burst  pressure  P4.  Also  shown  in 
this  figure  are  predictions  by  the  Stanford  code.  Agreement 
between  the  Stanford  model  and  experiment  is  reasonable. 


Fig.  15  Test  time  as  a  function  of  binary  gas  length  in  the  driver  section 
for  =  2.38.  As  Lbinarv  decreases,  the  test  time  increases.  Also, 

'^driven 

within  the  scatter  of  the  data,  this  relationship  is  approximately  linear. 
A  comparison  is  given  with  backfill  test  time  predictions  made  by  the 
Stanford  code  { =  1  .o)  as  well 

\  ^driven  / 


given  the  assumptions  employed  in  the  computer  algorithm. 
Notably,  the  model  correctly  predicts  the  approximate  rate 
of  test  time  increase  as  Lbinary  is  decreased.  The  underpre¬ 
diction  of  test  times  by  the  Stanford  code  can  be  attributed  in 
part  to  the  driver/driven  area  ratio  once  again;  at  this  temper¬ 
ature,  the  simulation  required  more  helium  in  the  binary  gas 
mixture,  thus  increasing  the  expansion  wave’s  velocity  and 
reducing  predicted  test  times.  Finally,  the  change  in  charac¬ 
ter  of  the  simulation  near  Lbinary  =  5  m  occurs  because  for 
shorter  Lbinary,  the  reflected  expansion  fan  wave  head  first 
encounters  the  reflected  shock,  whereas  for  longer  Lbinary, 
the  reflected  expansion  fan  wave  head  first  encounters  the 
backfill/binary  gas  interface. 


4.5  Variation  of  filling  time 

To  further  investigate  binary /backfill  gas  mixing,  the  flow  rate 
of  backfill  gas  was  varied;  this  resulted  in  different  filling 
times  for  the  backfill  gas.  Presumably,  as  the  backfill  gas 
flow  rate  is  increased  (i.e.,  shorter  backfill  time),  the  rate 
of  turbulent  mixing  between  the  binary  and  backfill  gases 
should  increase,  causing  the  binary /backfill  mixing  zone  to  be 
longer  and  affecting  the  resulting  test  time.  However,  in  our 
experiments  over  a  wide  range  of  backfilling  times,  the  test 
time  remained  constant.  Moreover,  the  non-ideal  pressure 
change  remained  approximately  constant  as  well.  This 
could  be  an  indication  that  the  length  of  the  mixing  zone  has 
little  effect  on  test  time,  or  that  other  factors  besides  backfill 
gas  flow  rate  are  more  important  in  determining  the  extent  of 
binary/backfill  gas  mixing. 
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4.6  Variation  of  backfill  gas  sound  speed 

To  investigate  the  effect  of  backfill  gas  sound  speed,  a  series 
of  experiments  were  conducted  using  different  backfill  gases. 
These  gases,  in  order  of  lowest  sound  speed  to  highest  sound 
speed  at  20  C,  were  100  %  CO2,  21  %  02/Ar,  100  %  Ar, 
21  %  O2/N2,  100  %  N2,  and  30  %  He/N2.  A  plot  relating 
test  time  to  backfill  gas  sound  speed  is  provided  in  Fig.  16. 
As  expected,  test  time  increases  as  sound  speed  decreases. 

4.7  Selection  of  backfill  gas 

The  selection  of  appropriate  backfill  gas  depends  on  a  num¬ 
ber  of  factors,  including  the  desired  test  time,  price  of  backfill 
gas,  component  vapor  pressure,  and  environmental  consid¬ 
erations  (several  potential  backfill  gases  and  their  properties 
are  provided  for  reference  in  Table  1).  First,  Fig.  16  can  be 
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Fig.  16  Test  time  as  a  function  of  backfill  gas  sound  speed  at  20  C. 
The  test  time  is  seen  to  increase  linearly  with  a  decrease  in  the  backfill 
gas  sound  speed 


used  as  a  selection  guide  for  determining  the  backfill  gas  nec¬ 
essary  to  obtain  the  desired  test  time,  holding  other  factors 
(e.g.,  Ldriver,  ^driven,  and  Lbinary)  Constant.  Second,  backfill 
gas  prices  are  somewhat  dependent  on  geographical  location. 
However,  in  general,  gases  obtained  from  the  atmosphere 
(e.g.,  N2  and  Ar),  natural  gas  (e.g.,  CH4  and  CsHg),  and 
combustion  products  (e.g.,  CO2)  have  reliably  low  prices. 
Third,  in  terms  of  gas  properties,  the  driver  section  burst 
pressure  (P4)  is  limited  by  the  vapor  pressure  of  the  backfill 
gas  [3].  For  example,  C3H8  and  CO2  have  vapor  pressures  of 
about  9  and  57  atm  at  20  C,  meaning  that  the  maximum  P4 
values  when  backfilling  with  these  species  are  9  and  57  atm, 
respectively.  Since  driver  burst  pressure  P4  generally  scales 
with  reflected  shock  pressure  L5  [22],  these  gases  may  be 
excluded  for  very  high-pressure  experiments.  Finally,  envi¬ 
ronmental  considerations  are  important  factors  in  choosing  a 
backfill  gas.  As  evident  from  Table  1,  sound  speed  generally 
decreases  with  gas  molecular  weight,  making  larger  mole¬ 
cules  attractive  candidates.  Several  examples  include  large 
hydrocarbons  (e.g.,  CsHg),  noble  gases  (e.g.,  Kr  and  Xe),  and 
insulators/refrigerants  (e.g.,  SFb  and  C4F10).  However,  these 
groups  possess  some  notable  drawbacks.  In  particular,  pres¬ 
surized  hydrocarbons  have  explosive  potential,  noble  gases 
are  often  prohibitively  expensive  unless  a  trace  gas  exhaust 
recovery  system  is  implemented,  and  insulators/refrigerants 
are  known  to  be  potent  greenhouse  gases.  Moreover,  some 
agents  such  as  CO  are  known  to  be  poisonous,  making  even 
small  leaks  very  dangerous.  As  such,  care  must  be  exercised 
in  selecting  the  most  appropriate  backfilling  gas  for  experi¬ 
mentation. 

5  Combination  of  techniques 

Figure  17  provides  a  demonstration  of  the  use  of  the  above 
techniques  in  concert.  This  experiment  employed  pure  argon 


0.00  0.01  0.02  0.03  0.04  0.05  0.06  0.07  0.08  0.09  0.10  0.11  seconds 


Fig.  17  Demonstration  of  combined  test  time  extension  techniques. 
This  shock  shows  102  ms  (0.102  s)  of  test  time  at  T5  =  524  K  and 
P5  =  1.6  atm;  the  dashed  line  shows  the  average  pressure  (1.695  atm). 
Conditions:  =  1-0,  Ldriven  =  9.73  m,  Ldriver  =  13.42  m, 


Lbinary  =  1.2  m.  Driver  100  %  N2  (at  this  temperature,  the  correct 
tailoring  requires  no  helium);  driven  100%  Ar;  backfill  100%  CO2.  No 
helium  interface  gas  was  needed.  A  series  of  driver  inserts  was  used  to 
mitigate  non-ideal  pressure  increases 
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in  the  9.73-m  driven  section.  The  13.42-m  driver  section  was 
filled  with  pure  nitrogen  gas  as  the  tailored  “binary”  gas  mix¬ 
ture  (0  %  He/100  %  N2)  and  pure  carbon  dioxide  as  the 
backfill  gas.  A  short  Lbinary  value  of  1.2  m  was  used.  More¬ 
over,  a  series  of  driver  inserts  were  employed  to  mitigate 
non-ideal  pressure  changes.  These  included  a  conical  insert 
similar  to  that  of  Ref.  [6]  to  reduce  the  driver  section  diame¬ 
ter,  such  that  =  1.0,  allowing  the  reflected  shock  wave 
to  pass  back  into  the  driver  section  with  no  immediate  area 
change.  (For  >  1,  an  expansion  wave  is  generated  and 
propagates  back  to  the  test  location  when  the  reflected  shock 
passes  back  into  the  driver  section;  likewise,  for  <  1 ,  a 

■^driven 

shock  wave  is  generated  and  propagates  back  to  the  test  sec¬ 
tion).  In  this  experiment,  the  resulting  post-reflected  shock 
temperature  was  sufficiently  low  that  no  helium  interface  gas 
was  required  as  discussed  in  Ref.  [4].  The  resulting  pressure 
trace  shows  a  test  time  of  102  ms.  To  our  knowledge,  this  is 
the  longest  constant-pressure  test  time  to  have  been  achieved 
in  a  reflected  shock  experiment. 

Figure  18  shows  the  Stanford  code’s  estimate  as  to  the 
test  time  obtainable  as  a  function  of  T3  with  and  with¬ 
out  backfilling,  for  L driven  =  9.73  m,  L driver  =  13.42  m. 
Lbinary  =  1.2  m,  tailored  He/N2  driver  gas,  argon  driven  gas, 
and  carbon  dioxide  backfill  gas.  Importantly,  the  introduction 
of  staged  driver  gas  filling  can  substantially  increase  the  avail¬ 
able  test  time  at  all  temperatures.  Relative  to  the  experiment 
shown  in  Fig.  17,  the  model  underpredicts  the  measured  test 
time.  Again,  this  is  likely  due  to  the  idealized  assumptions 
implicit  in  the  computer  code,  such  as  zero  backfill/binary 
gas  mixing,  but  could  also  be  the  result  of  uncertainty  in  the 
experimental  value  of  Lbinary  =  1.2  m,  since  this  value  does 
not  account  for  the  volume  occupied  by  the  inserts  employed. 
Note  that  the  code  predicts  that  backfilling  allows  for  greater 


Fig.  18  Benefit  of  backfill  technique  on  test  time  versus  T5 .  Lines  show 
Stanford  code  predictions  and  blue  square  shows  experimental  measure¬ 
ment.  Conditions:  =  1.0,  Ldriven  =  9.73  m,  Ldriver  =  13.42  m, 

'/I  driven 

I-binary  =  1.2  m.  Driver  He/N2  (tailored);  driven  100  %  Ar;  backfill 
100  %  CO2.  No  helium  interface  gas 


increases  in  test  time  at  higher  temperatures  as  compared  to 
low  temperatures.  This  reflects  the  fact  that  at  high  tempera¬ 
tures,  a  majority  of  the  tailored  driver  gas  consists  of  helium, 
implying  that  replacing  this  gas  with  a  low- sound- speed  alter¬ 
native  can  yield  great  performance  enhancements. 

6  Sample  ignition  delay  time  measurement 

We  turn  now  to  an  example  application  of  the  long  constant- 
pressure  test  time  concept  by  examining  a  problem  of 
practical  interest,  namely  the  measurement  of  ignition  delay 
time.  A  sample  ignition  delay  time  measurement  for  an  ^- 
heptane/21  %  oxygen/argon  (0  =  0.5)  mixture  at  4.6  atm 
and  648  K  is  provided  in  Fig.  19.  Details  of  the  diagnos¬ 
tics  employed  in  this  experiment  are  available  elsewhere 
[15],  hence  only  a  concise  summary  will  be  provided  here. 
At  these  conditions,  ^-heptane  exhibits  two-stage  ignition. 
The  first- stage  ignition  was  monitored  by  observing  the  fuel 
decay  according  to  the  Beer-Lambert  relation  using  a  contin¬ 
uous  wave  helium  neon  (HeNe)  laser  at  3.39  |Jim  that  probes 
Region  5  at  a  location  2  cm  from  the  endwall;  the  decline 
in  absorbance  near  50  ms  clearly  marks  this  event.  Second- 
stage  ignition  is  monitored  by  emission  from  excited  OH 
(OH*)  radicals  near  306  nm  through  a  port  in  the  shock  tube 
endwall;  the  rise  in  the  emission  signal  near  65  ms  indicates 
this  occasion  and  is  corroborated  by  the  rise  in  pressure  and 
the  final  drop  in  HeNe  absorbance.  The  second- stage  ignition 
delay  time  shown  is  >65  ms,  making  it,  to  our  knowledge, 
one  of  the  longest  ignition  delay  times  ever  measured  in  a 


Fig.  19  Sample  ignition  delay  time  measurement.  The  first  drop  in 
fuel  absorbance  indicates  first-stage  ignition;  the  rise  in  OH*  emission, 
rise  in  pressure,  and  further  drop  in  fuel  absorbance  indicate  second- 
stage  ignition.  Conditions:  =  1.0,  Ldriven  =  9.73  m,  Ldriver  = 

12.61  m,  Lbinary  =  4.2  m.  Driver  34  %  He/N2;  driven  n-heptane/21  % 
02/Ar;  backfill  100  %  N2.  A  driver  insert  was  employed  to  eliminate 
non-ideal  pressure  increases,  and  no  helium  interface  gas  was  used. 
rign,i  =  48.0  ms;  Tign,2  =  65.8  ms 
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shock  tube.  This  experiment  would  not  have  been  possible 
without  the  test  time  extension  techniques  discussed  in  this 
work  and  the  results  provide  a  good  example  of  the  nega¬ 
tive  temperature  coefficient  (NTC)  region  behavior  exhibited 
by  many  large  alkane  fuels  at  practical  engine  conditions.  It 
is  worth  noting  that  these  experiments  provide  data  similar 
to  that  obtainable  in  a  rapid  compression  machine  (RCM), 
but  without  ambiguity  associated  with  the  finite  compression 
time  and  gradual  pressure  decays  in  such  facilities. 


7  Heat  transfer 


As  mentioned  in  Sect.  1.2,  a  numerical  analysis  by  Frazier 
et  al.  [18]  concluded  that  wall  heat  transfer  does  not  reduce 
area- averaged  temperature  at  long  times,  for  shock  tubes  with 
8  cm  or  larger  internal  diameters.  To  examine  this  statement, 
we  measured  gas  pressure  and  temperature  in  situ  for  several 
carbon  dioxide- seeded  experiments  (4  %  C02/Ar)  over  long 
test  durations.  Like  the  ignition  measurement,  details  of  the 
CO2 -based  spectroscopic  temperature  diagnostic  are  avail¬ 
able  elsewhere  [15,30],  so  only  a  brief  overview  is  provided 
here.  The  temperature  diagnostic  was  based  on  a  two-color 
absorption  technique,  which  exploited  the  steep  inverse  tem¬ 
perature  dependence  between  CO2  transitions  in  different 
infrared  absorption  bands  to  achieve  high  sensitivity.  Spec¬ 
troscopic  measurement  uncertainty  is  estimated  to  be  <2  % 
for  the  reflected  shock  conditions  in  these  experiments. 

Figure  20  depicts  one  of  these  shocks,  taken  at  635  K  and 
5.1  atm.  During  this  experiment,  a  very  small  residual  non¬ 
ideal  pressure  change  of  =0.22  %/ms  was  allowed  to 
examine  the  isentropic  relationship  between  gas  temperature 
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Fig.  20  Sample  CO2 -based  spectroscopic  temperature  measurement 
and  comparison  with  isentropically  computed  temperature  (1).  Condi¬ 
tions:  ^47^^  —  2.38,  L^j-iven  =  9.73  m,  L^nver  —  9.57  m,  L^inary  — 

^driven 

6.2  m.  Driver  6.3  %  He/N2;  driven  4  %  C02/Ar  (ys  =  1.6128);  back¬ 
fill  100  %  N2.  A  driver  insert  was  employed  to  control  the  non-ideal 
pressure  increase,  and  no  helium  interface  gas  was  used 


and  pressure  ( 1 );  near  the  end  of  the  test  time  this  amounted  to 
an  1 1  %  overall  increase  in  pressure.  Several  important  notes 
can  be  made  about  this  graphic.  First,  the  initial  post-reflected 
shock  temperature  computed  based  on  incident  shock  speed 
matches  the  spectroscopically  obtained  value  to  within  0.5  %, 
which  is  the  approximate  uncertainty  in  the  measured  tem¬ 
perature;  this  underscores  the  reliability  of  the  velocity-based 
calculation  technique  used  widely  in  shock  tube  studies. 
Moreover,  the  uncertainty  in  shock  speed-based  temperature 
for  this  experiment  is  conservatively  estimated  to  be  only 
about  0.3  %  by  accounting  for  errors  in  measured  incident 
shock  speed  (0.14  %  in  this  trial),  initial  (pre-shock)  tem¬ 
perature  and  pressure,  and  thermodynamic  properties  (heat 
capacity  and  enthalpy).  At  this  stage  of  development  of  the 
laser-based  diagnostic,  there  is  slightly  less  uncertainty  in 
the  shock  speed-based  temperature  value  than  that  mea¬ 
sured  spectroscopically,  but  the  important  message  is  that 
these  values  agree  within  their  (small)  experimental  uncer¬ 
tainties.  Second,  the  subsequent  spectroscopically-measured 
temperature  change  matches  that  which  was  calculated  from 
the  change  in  pressure  [assuming  isentropic  compression 
according  to  (1)]  very  closely.  This  gives  confidence  for 
use  of  the  isentropic  relationship  to  estimate  changes  in 
post-reflected  shock  temperature  based  on  observed  pressure 
profiles.  Hence,  the  conclusion  of  Frazier  et  al.  that  wall  heat 
transfer  is  not  a  significant  concern  in  large-diameter  shock 
tubes,  like  that  employed  in  this  work,  at  conditions  similar  to 
those  of  this  study,  seems  well  founded  [18].  Were  heat  trans¬ 
fer  a  problem,  we  might  expect  a  significant  temperature  and 
pressure  drop  over  the  course  of  the  experiment.  Note  that 
the  spectroscopic  diagnostic  yields  a  path-average  tempera¬ 
ture  value;  due  to  conduction  heat  transfer,  the  temperature 
of  the  gas  in  the  thin  boundary  layer  at  the  edges  of  the  shock 
tube  will  be  lower  than  the  bulk  temperature.  Further  work 
is  needed  to  investigate  the  role  of  heat  transfer  over  a  wider 
range  of  conditions,  including  higher  temperatures. 

Finally,  it  should  be  noted  that  under  some  conditions, 
shock  bifurcations  can  cause  additional  effects.  Shock  bifur¬ 
cations  occur  as  the  reflected  shock  interacts  with  the  gas 
boundary  layer  following  the  incident  shock.  While  a  com¬ 
plete  treatment  of  the  effects  of  shock  bifurcations  and 
boundary  layers  is  beyond  the  scope  of  this  work,  it  should 
be  noted  that  these  issues  are  less  significant  when  experi¬ 
menting  with  large-bore  shock  tubes,  monatomic  gases,  and 
low  pressures  [31-33].  Additionally,  the  effects  of  the  shock 
bifurcation  on  the  post-reflected  shock  gases  are  most  impor¬ 
tant  immediately  following  the  reflected  shock  wave  and  not 
at  long  times  (though  boundary  layer  mixing  that  affects  the 
shock  tube  temperature  field  can  still  occur  at  long  timescales 
like  those  demonstrated  in  this  work).  Nevertheless,  the  spec¬ 
troscopic  temperature  measurements  obtained  in  this  work 
lend  confidence  to  the  use  of  shock  tubes  at  long  test  times 
with  well-defined  thermodynamic  conditions. 
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8  Concluding  remarks 

We  have  demonstrated  the  use  of  several  shock  tube  tech¬ 
niques  for  achieving  extended  reflected  shock  test  times  with 
uniform  pressure,  including  driver  inserts,  driver  gas  tailor¬ 
ing,  the  use  of  helium  interface  gas  in  front  of  the  diaphragm, 
driver  extensions,  and  staged  driver  gas  filling.  The  latter 
method,  which  involves  first  filling  the  driver  section  with 
a  tailored  binary  gas  mixture,  and  secondly  backfilling  a 
low-speed-of-sound  gas  near  the  driver  section  end  cap,  can 
increase  test  times  by  approximately  20  %  while  reduc¬ 
ing  necessary  driver  helium  consumption  by  up  to  about 
85  %.  By  combining  the  above  techniques  together,  we  have 
achieved  an  experiment  with  102  ms  of  test  time  at  524  K  and 
1.6  atm.  We  have  demonstrated  the  value  of  these  strategies 
by  measuring  long  ^ -heptane/21  %  02/Ar  (0  =  0.5)  igni¬ 
tion  delay  times  at  pressures  near  4.5  atm  and  temperatures 
near  650  K.  Also,  we  have  employed  a  novel  spectroscopic 
CO2 -based  temperature  diagnostic  to  demonstrate  the  via¬ 
bility  of  using  isentropic  relations  to  predict  changes  in  test 
temperature  based  on  the  measured  pressure  profile  for  each 
shock,  and  hence  have  shown  that  temperature  remains  con¬ 
stant  throughout  the  entire  test  duration  in  long  shocks  when 
pressure  is  held  constant.  The  techniques  introduced  herein 
will  allow  use  of  shock  tubes  to  explore  low-temperature, 
moderate-pressure  chemistry  which  is  of  practical  impor¬ 
tance  to  a  variety  of  combustion  devices  in  use  today. 
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Abstract 

Ignition  delay  times  of  normal  heptane  have  been  measured  at  temperatures  ranging  from  65 1  to  823  K 
and  at  pressures  between  6.1  and  7.4  atm  at  an  equivalence  ratio  of  0.75  in  15%02/5%C02/Ar  and  in 
15%02/Ar  mixtures  behind  reflected  shock  waves  in  a  shock  tube.  Time-history  measurements  of  fuel, 
OH,  aldehydes  (mostly  CH2O),  CO2,  H2O,  and  temperature  were  also  measured  under  these  conditions. 
These  time-histories  provide  critically  needed  kinetic  targets  to  test  and  refine  large  reaction  mechanisms. 
Measurements  were  acquired  using  a  novel  constrained  reaction  volume  approach,  wherein  a  sliding  gate 
valve  confined  the  reactant  mixture  to  a  region  near  the  endwall  of  the  shock  tube.  A  staged-driver  gas  fill¬ 
ing  strategy,  combined  with  driver  section  extensions,  driver  inserts,  and  driver  gas  tailoring,  was  used  to 
obtain  constant-pressure  test  times  of  up  to  55  ms,  allowing  observations  of  the  chemistry  in  the  Negative 
Temperature  Coefficient  (NTC)  region.  Experiments  with  conventional  shock  tube  filling  were  also  per¬ 
formed,  showing  similar  overall  ignition  behavior.  Comparisons  between  current  data  and  simulations 
using  the  Mehl  et  al.  /i-Heptane  mechanism  (2011)  are  provided,  revealing  that  the  mechanism  generally 
under-predicts  first-stage  ignition  delay  times  in  the  NTC  region,  and  that  at  low  temperatures  it  over¬ 
predicts  the  extent  of  fuel  decomposition  during  first  stage  ignition. 

©2014  The  Combustion  Institute.  Published  by  Elsevier  Inc.  All  rights  reserved. 
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1.  Introduction 

Recent  efforts  to  improve  internal  combustion 
engine  efficiency  have  resulted  in  the  need  for 
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better  kinetic  mechanisms  that  describe  fuel  oxi¬ 
dation  chemistry  at  practical  conditions.  One 
important  fuel  that  has  been  the  subject  of  numer¬ 
ous  experimental  (e.g.,  [1-5])  and  modeling  (e.g., 
[6,7])  studies  is  normal  heptane,  as  it  is  a  primary 
reference  fuel. 

To  date,  most  shock  tube  ignition  delay  time 
studies  of  22-heptane  at  Negative  Temperature 
Coefficient  (NTC)  region  temperatures  (between 
600  and  900  K)  (e.g.,  [1-3])  have  been  performed 
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at  pressures  above  10  atm,  since  available  test 
times  for  shock  tubes  have  historically  been 
restricted  to  about  10  ms  and  ignition  delay  times 
are  longer  than  this  at  lower  pressures.  As  a  result, 
leading  kinetic  mechanisms  (e.g.,  that  of  Lawrence 
Livermore  National  Laboratory,  LLNL  [7])  have 
not  been  validated  at  NTC  region  temperatures 
and  moderate  pressures.  Moreover,  long  ignition 
delay  times  experimentally  measured  in  shock 
tubes  have  been  subject  to  pre-ignition  phenom¬ 
ena  occurring  away  from  the  test  location,  which 
have  been  known  to  cause  an  increase  in  pressure 
and  temperature  and  ultimately  decrease  mea¬ 
sured  ignition  delay  times  [2].  This  creates  chal¬ 
lenges  in  simulating  measured  ignition  delay 
times,  as  conditions  during  each  experiment  can¬ 
not  be  properly  classified  as  possessing  constant 
internal  energy  and  constant  volume  (U-V)  [8,9] 
as  has  been  commonly  assumed. 

The  need  for  experimental  data  that  are  easy  to 
simulate  at  engine-relevant  (moderate-pressure/ 
low-temperature)  conditions  is  evident.  Thus  we 
have  designed  a  new  shock  tube  that  is  able  to 
achieve  near-constant-pressure  test  times  of  up 
to  55  ms  using  a  family  of  techniques  including 
driver  extensions  [10],  driver  inserts  [11],  driver 
gas  tailoring  [12],  and  staged-driver  gas  filling 
[13].  Moreover,  this  shock  tube  employs  a  sliding 
gate  valve  to  confine  the  reactant  mixture  to  a 
region  near  the  endwall,  minimizing  total  energy 
release  (and  hence  pressure  change)  and  eliminat¬ 
ing  the  effects  of  chemistry  occurring  far  from  the 
endwall;  we  call  this  implementation  a  constrained 
reaction  volume  (CRV)  shock  tube  [8,9,14].  Using 
this  device,  we  have  measured  ^-heptane  ignition 
delay  times,  mole  fraction  time-histories  (fuel, 
OH,  aldehydes,  CO2,  and  H2O),  and  temperature 
histories  at  NTC  region  conditions. 


2.  Experimental  setup 

2.1.  Shock  tube  facility 

The  new  CRV  shock  tube  consists  of  a  9.57-m 
driver  section  filled  with  a  high-pressure  tailored 
mixture  of  helium  and  nitrogen  (see  [12]),  sepa¬ 
rated  by  a  polycarbonate  diaphragm  from  a  low- 
pressure  mixture  called  the  buffer  gas  (CO2/N2/ 
Ar)  in  a  9.32-m  driven  section  (internal  diameter 
of  11.53  cm).  Positioned  at  the  end  of  the  buffer 
gas  section  is  the  sliding  gate  valve,  which  confines 
the  test  gas  mixture  to  a  39-cm  long  test  section 
adjacent  to  the  endwall.  The  buffer  gas  is  specially 
formulated  so  that  its  molecular  weight  and  spe¬ 
cific  heat  ratio  7  =  ^  match  that  of  the  test  gas 
at  the  initial  temperature  of  the  experiment.  Addi¬ 
tionally,  the  driver  section  is  equipped  with  two 
filling  ports  to  allow  for  staged-driver  gas  filling; 
one  port  is  used  to  initially  fill  the  tailored 
helium/nitrogen  mixture,  and  the  second,  located 


far  from  the  diaphragm,  allows  a  low-sound-speed 
gas  (in  these  experiments,  nitrogen)  to  be  filled 
subsequently.  Since  the  secondary  port  is  located 
near  the  end  of  the  long  driver  section,  the  tailored 
helium/nitrogen  mixture  immediately  adjacent 
to  the  diaphragm  remains  undiluted  prior  to 
experimentation. 

A  shock  experiment  commences  when  the  final 
nitrogen  is  added  through  the  secondary  driver 
port,  causing  the  diaphragm  to  burst.  This  rapid 
bursting  results  in  the  formation  of  an  incident 
shock  wave,  which  moves  toward  the  test  section, 
and  an  expansion  wave,  which  moves  back  into 
the  driver  section.  The  velocity  of  the  shock  wave 
is  measured  using  a  series  of  seven  PCB  model 
113A26  transducers.  This  velocity  is  extrapolated 
to  the  shock  tube  endwall  with  an  accuracy  of 
±0.2%  and  is  used  to  calculate  post-shock  temper¬ 
ature  and  pressure  to  within  ±0.6%  and  ±1.5%, 
respectively,  using  standard  gasdynamic  relations 
and  thermodynamic  properties.  Pressure  is  mea¬ 
sured  in  situ  using  a  Kistler  model  601 B1  trans¬ 
ducer  located  in  the  test  section,  confirming 
computed  conditions.  A  few  seconds  prior  to  dia¬ 
phragm  rupture,  the  sliding  gate  valve  is  opened; 
in  situ  spectroscopic  measurements  of  fuel  concen¬ 
tration  have  shown  that  mixing  between  the  buffer 
and  test  gases  during  this  short  interval  is  negligi¬ 
ble.  Since  the  buffer  and  test  gases  have  nearly 
identical  molecular  weights  and  specific  heat 
ratios,  the  velocity  and  attenuation  rate  of  the 
shock  wave  are  constant  across  the  buffer-test 
gas  interface.  Moreover,  this  formulation  prevents 
reflected  shock/expansion  waves  at  this  interface. 
As  the  shock  wave  propagates  into  the  test  gas, 
it  heats  and  compresses  this  gas  from  its  original 
volume,  spanning  a  distance  Li  =  39  cm  between 
the  endwall  and  the  gate  valve,  to  a  smaller  vol¬ 
ume  close  to  the  endwall.  Ultimately,  the  shock 
wave  reflects  from  the  endwall  of  the  test  section 
and  moves  back  through  the  test  gas  toward  the 
diaphragm,  further  heating  and  compressing  the 
test  gas.  This  results  in  the  test  gas  being  stagnated 
and  confined  in  a  volume  with  a  nominal  length 
±5  as  measured  from  the  endwall,  enabling  spec¬ 
troscopic  and  pressure  diagnostics  to  probe  the 
reactive  test  gas  at  a  location  2.0  cm  from  the 
endwall.  Typical  values  of  L5  in  these  CRV  exper¬ 
iments  were  between  5  and  10  cm,  whereas  values 
of  ±5  for  conventional  (entire  shock  tube)  test  gas 
filling  were  between  130  and  240  cm.  Additional 
details  concerning  the  CRV  shock  tube  technique 
can  be  found  in  Ref.  [14]. 

Meanwhile,  the  expansion  wave  in  the  driver 
section  proceeds  through  the  tailored  helium/ 
nitrogen  mixture  and  then  propagates  into  the 
back-filled  pure  nitrogen  at  the  end  of  the  driver 
section.  Since  this  expansion  wave  travels  at  the 
local  speed  of  sound,  its  velocity  is  greatly  reduced 
in  the  pure  nitrogen  compared  to  that  in  the 
tailored  helium/nitrogen  mixture.  Eventually  this 
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expansion  wave  reflects  from  the  endwall  of  the 
driver  section  and  proceeds  back  toward  the  test 
section;  its  arrival  at  the  test  section  ends  the 
experiment.  The  use  of  staged-driver  filling 
substantially  increases  test  time  by  delaying  this 
arrival.  Additional  details  concerning  the  staged- 
driver  gas  filling  technique  can  be  found  in  Ref. 

Mixtures  were  prepared  manometrically  in 
stainless  steel  mixing  tanks  pumped  to  less  than 
0.10  Torr  prior  to  gas  mixing.  The  shock  tube 
was  turbo-mechanically  pumped  to  better  than 
10“^  Torr  prior  to  most  experiments,  with  negligi¬ 
ble  subsequent  leak-plus-outgassing.  Normal  hep¬ 
tane  was  obtained  from  Sigma  Aldrich  (99.3%) 
and  was  degassed  prior  to  use;  research-grade 
oxygen,  carbon  dioxide,  nitrogen,  and  argon  were 
obtained  from  Praxair,  Inc. 

2.2.  Laser  absorption  measurements 

Details  of  the  laser  absorption  techniques 
applied  in  this  study  are  available  in  the  papers 
referenced  in  this  section,  hence  only  a  brief  treat¬ 
ment  will  be  given  here.  Laser  diagnostics  are  able 
to  make  quantitative  measurements  of  gaseous 
species  because  the  ratio  of  monochromatic  light 
intensity  at  frequency  v  incident  upon  an  absorb¬ 
ing  gas  mixture  4  to  that  exiting  the  mixture  /  is 
related  to  the  number  density  of  absorbing  species 
n  according  to  the  Beer-Lambert  relation. 


In  this  equation,  a  is  known  as  the  absorbance,  Oy 
is  the  propensity  for  light  absorption  by  the  mol¬ 
ecule  at  frequency  v  (called  the  absorption  cross 
section)  and  L  is  the  path  length  along  which 
the  absorption  occurs. 

2.2.1.  Laser  absorption  of  fuel 
Measurements  of  fuel  concentration  were 

made  using  the  3.39  pm  output  of  a  continuous 
wave  (CW)  Helium-Neon  (HeNe)  laser.  Cross 
section  values  were  obtained  from  [15-17]  and 
were  confirmed  using  shock  tube  experiments. 
Before  each  shock  experiment,  the  mole  fraction 
of  fuel  was  measured  using  this  spectroscopic 
technique;  measured  values  were  typically  3% 
lower  than  manometric  values,  attributed  to  fuel 
loss  to  the  walls  of  the  unheated  shock  tube.  Fuel 
mole  fraction  values  were  also  measured  between 
the  incident  and  reflected  shocks  on  each  experi¬ 
ment,  confirming  that  mixing  of  buffer  and  test 
gases  had  not  diluted  the  reactive  mixture. 

2.2.2.  Laser  absorption  of  OH  and  aldehydes 
OH  concentration  was  measured  using  the  out¬ 
put  of  a  narrow-linewidth  CW  laser  at  306.7  nm, 
the  peak  of  the  well-characterized  Ri(5)  absorp¬ 


tion  line  in  the  OH  A-X  (0,0)  band  [18,19].  The 
laser  light  was  obtained  by  intracavity  fre¬ 
quency-doubling  of  613.4  nm  light  generated  by 
a  Spectra  Physics  380  ring  dye  laser  that  was 
pumped  by  a  5  W  Coherent  Verdi  at  532  nm. 
Common  mode  rejection  was  used  to  reduce  laser 
intensity  noise,  resulting  in  a  minimum  detectivity 
of  5  ppm  at  the  conditions  of  this  work,  with  an 
estimated  uncertainty  of  ±10%.  Measurements 
were  also  performed  with  the  laser  tuned  off  the 
OH  line  to  observe  other  absorbing  species, 
enabling  use  of  a  simple  subtraction  scheme  to 
cancel  spectrally  broad  interference.  At  the  condi¬ 
tions  of  these  experiments,  interference  is  largely 
due  to  molecules  possessing  a  C=0  bond,  i.e., 
aldehydes,  ketenes,  and  formic  acid.  Kinetic  simu¬ 
lations  performed  using  the  LLNL  mechanism  [7] 
revealed  that  CH2O  is  the  dominant  aldehyde  pro¬ 
duced  during  first  stage  ignition  of  ^-heptane. 
Therefore,  the  cross  section  of  CH2O  [20]  was 
assumed  to  represent  that  of  aldehydes/ketenes/ 
formic  acid,  and  time-histories  of  these  species’ 
total  concentration  have  thus  been  inferred  and 
reported  here.  Finally,  emission  by  OH  at 
306  nm  in  the  laser  absorption  scheme,  which  used 
a  series  of  irises  to  limit  radiation  reaching  the 
detector,  was  found  to  be  negligible. 

2.2.3.  Laser  absorption  of  CO 2 

Carbon  dioxide  mole  fraction  measurements 
were  attained  by  a  fixed-wavelength  direct  absorp¬ 
tion  method  near  4.2  pm  using  a  CW  external- 
cavity  quantum  cascade  laser  [21].  The  laser  was 
tuned  to  the  R(66)  transition  in  the  fundamental 
band  near  2387.24  cm“\  which  offers  minimal 
absorption  temperature  sensitivity  over  the  tem¬ 
perature  range  700-1000  K.  The  measurement 
uncertainty  is  <10%,  and  the  minimum  detectivity 
is  approximately  5  ppm. 

2.2.4.  Laser  absorption  of  H2O  and  temperature 
measurement 

Water  mole  fraction  and  temperature  measure¬ 
ments  were  acquired  along  a  single  line-of-sight 
using  two  frequency-multiplexed  tunable  diode 
lasers  centered  upon  H2O  absorption  transitions 
near  3920.09  and  4029.52  cm“^  The  wavelength 
of  each  laser  was  scanned  across  the  transition 
linecenter  at  20  kHz  and  modulated  at  90  and 
112  kHz,  respectively,  to  monitor  first-harmonic- 
normalized  wavelength-modulation  spectroscopy 
with  second-harmonic  detection  (WMS-2//1/)  sig¬ 
nals  at  linecenter.  Details  regarding  the  experi¬ 
mental  setup  and  select  spectroscopic  parameters 
are  given  in  [22,23]. 

2.3.  Light  emission 

The  emission  from  kinetically  excited  OH  rad¬ 
icals  (OH*)  at  306  ±75  nm  was  also  measured 
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during  ignition  using  a  silicon  photo-detector  and 
a  UG-5  Shott  glass  band-pass  filter.  This  diagnos¬ 
tic  was  positioned  at  both  a  window  located  2  cm 
from  the  endwall  and  at  a  port  in  the  endwall; 
ignition  delay  time  measurements  at  the  two  loca¬ 
tions  generally  agreed  to  within  10%. 

2.4.  Ignition  measurement 

Ignition  delays  were  measured  relative  to  the 
time  of  arrival  of  the  reflected  shock,  taken  here 
as  the  time  of  the  half-rise  in  pressure  caused  by 
the  reflected  shock.  First  and  second  stage  ignition 
events  in  the  NTC  region  were  measured  by 
observing  the  fuel  decay  process  by  way  of  the 
3.39  |im  HeNe  diagnostic.  The  time  of  first  stage 
ignition  was  taken  to  be  the  instant  of  maximum 
slope  in  the  initial  decline  of  the  fuel  absorbance, 
and  the  time  of  second  stage  ignition  was  taken  to 
be  the  instant  when  the  absorbance  declined  to  5% 
of  its  initial  post-reflected-shock  value.  At  high 
temperatures,  the  maximum  OH*  emission  slope 
was  extrapolated  to  the  baseline  (zero)  value  to 
determine  the  ignition  delay  time. 

A  sample  ignition  measurement,  including  a 
comparison  between  CRV  and  conventional 
shock  tube-filling  experimentation,  is  shown  in 
Fig.  1.  Observe  that  the  decreasing  HeNe  absor¬ 
bance  provides  a  clear  indication  of  the  two  stages 
of  ignition.  Moreover,  comparing  the  reactive 
pressure  traces  to  the  non-reactive  trace  near  these 
ignition  events  serves  to  illustrate  the  differences 
between  CRV  and  conventional  experimentation. 
Near  the  first  ignition  event,  the  conventional 
shock  exhibits  a  clearly  visible  and  sustained  pres¬ 
sure  increase,  whereas  the  CRV  experiment  shows 
virtually  no  rise  relative  to  the  nonreactive  pres- 


Fig.  1.  Comparison  between  CRV  and  conventional 
shock  measurements.  Reflected  shock  initial  conditions: 

^  0.75  («-heptane/l  5%02/5%C02/Ar),  6.6  atm, 

~703  K,  Ls  =  9  cm.  Ignition  delay  times  (indicated  by 
open  and  solid  circles):  ti,crv  =  8.2  ms,  ti,crv  =  26.8  - 

conventional  8.8  mS,  T2,conventional  22.7  mS. 


sure.  Furthermore,  at  the  second  ignition  event, 
the  conventional  shock  shows  an  exponential 
pressure  rise  to  a  level  roughly  four  times  that  of 
the  non-reactive  experiment,  whereas  the  CRV 
experiment  shows  only  a  small  pressure  increase 
that  quickly  relaxes  to  the  non-reactive  pressure. 
Importantly,  since  the  CRV  shock  remains  at 
nearly  constant  pressure  up  to  and  even  through¬ 
out  its  ignition  process,  this  experiment  can  be 
effectively  modeled  as  a  constant  enthalpy-pres¬ 
sure  (H-P)  reactor.  In  contrast,  the  conven¬ 
tional-fill  experiment  cannot  be  correctly 
modeled  in  a  simple  way,  e.g.  as  a  zero-dimen¬ 
sional  reactor  at  constant  energy-volume  (U-V), 
though  such  an  approach  may  be  adequate  in 
many  cases  for  simulating  first  stage  ignition,  or 
in  strong  ignition  cases  with  no  pre-ignition 
energy  change.  The  slight  difference  in  first- stage 
ignition  time,  between  conventional  and  CRV 
experiments,  is  likely  due  to  small  mismatches  in 
equivalence  ratio  and  initial  reaction  temperature. 


3.  Results 

3.1.  Ignition  delay  times 

Ignition  delay  times  for  ^-heptane  were  mea¬ 
sured  using  the  CRV  strategy  at  an  equivalence 
ratio  of  0.75  in  15%02/5%C02/Ar  and  in 
15%02/Ar  mixtures  at  6.5  atm  and  temperatures 
between  651  and  823  K.  (CO2  was  included  in 
mixtures  in  support  of  a  new  C02-absorption 
based  temperature  diagnostic,  not  reported  here. 
Additionally,  CO2  was  not  included  in  the  deter¬ 
mination  of  equivalence  ratio,  and  simulations 
[7]  revealed  that  the  presence  of  CO2  had  negligi¬ 
ble  effects  on  ignition  kinetics).  Measurements 
were  also  conducted  using  the  conventional  filling 
strategy  throughout  this  range  for  purposes  of 
comparison.  In  addition  to  these  low  temperature 
tests,  experiments  were  conducted  using  both  the 
CRV  and  conventional  strategies  at  higher  tem¬ 
peratures  and  a  pressure  of  7  atm.  These  ignition 
delay  data  are  summarized  in  Fig.  2.  Small  scaling 
corrections,  obtained  from  the  LLNL  mechanism 
[7]  by  perturbing  pressure  and  equivalence  ratio  at 
several  temperatures  and  observing  the  effect  on  ti 
and  T2,  to  ignition  delay  times  have  been  applied 
to  normalize  individual  points  to  a  common  P 
and  (f)  (these  scaling  factors  are  temperature- 
dependent).  Furthermore,  ignition  delay  time  pre¬ 
dictions  made  using  the  LLNL  [7]  mechanism 
have  been  included  in  this  figure. 

Several  observations  can  be  made  about  Fig.  2. 
First,  compare  the  constant  U-V  and  constant  H- 
P  modeling  methods  (henceforth  referred  to  as 
simply  “UV”  and  “HP”,  respectively).  Observe 
that  at  high  temperatures,  the  UV  and  HP  simula¬ 
tions  yield  nearly  identical  results.  However,  in  the 
NTC  region,  the  UV  simulation  yields  strikingly 
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Temperature  (1000/K) 


Fig.  2.  w-Heptane  ignition  delay  time  values  (^-heptane/ 
15%02/5%C02/Ar,  0  =  0.75,  6.5  atm).  Solid  lines: 
constant  H-P  modeling;  dashed  lines:  constant  U-V 
modeling.  NTC  region  data  are  presented  at  the  initial 
post-reflected  shock  temperature,  prior  to  the  tempera¬ 
ture  increase  associated  with  flrst  stage  ignition.  Error 
bars  are  not  visible  at  the  resolution  used  in  this  plot. 


lower  second  stage  (i.e.,  complete)  ignition  delay 
times  than  the  HP  simulation,  while  the  two  pre¬ 
dict  identical  first  stage  times.  This  paradox  is 
the  result  of  an  ignition-accelerating  pressure 
increase  which  occurs  in  UV  simulations  during 
first  stage  ignition  but  which  does  not  occur  in 
HP  simulations.  Since  the  increase  occurs  during 
the  first  stage  ignition  but  not  before,  differences 
between  UV  and  HP  trials  exhibit  themselves  in 
the  second  stage  ignition  delay  times  rather  than 
in  the  first. 

Next,  compare  the  CRV  and  conventional 
shock  tube  experiments.  Notice  that  at  high  tem¬ 
peratures,  CRV  and  conventional  data  agree 
closely  with  one  another  (as  well  as  with  simula¬ 
tions);  this  similarity  was  expected,  as  nonideal 
effects  in  shock  tubes  typically  become  important 
at  long  times,  i.e.  at  lower  temperatures  [10].  Also 
in  agreement  are  the  CRV  and  conventional  first 
stage  ignition  delay  times  in  the  NTC  (between 
825  and  650  K),  consistent  with  the  fact  that 
CRV  and  conventional  pressure  traces  looked 
identical  up  to  the  point  of  first  stage  ignition. 
Furthermore,  at  temperatures  below  about 
725  K,  the  CRV  and  conventional  techniques 
yield  similar  second  stage  ignition  delay  times  with 
low  scatter,  in  contrast  with  expectation  based  on 
simulations.  However,  between  about  825  and 
725  K,  the  second  stage  ignition  data  show  signif¬ 
icantly  more  scatter;  one  possible  explanation  is 
that  ignition  phenomena  at  these  temperatures 
are  more  sensitive  to  the  small  deviations  in  pres¬ 
sure  present  in  these  experiments  relative  to  the 
nominal  value  of  6.5  atm.  The  LENT  mechanism 
[7]  supports  this  hypothesis;  it  gives  a  pressure 
scaling  of  at  800  K,  whereas  at  675  K  it  is 


P“^'^.  In  general,  it  appears  that  conventional 
shocks  performed  between  825  and  725  K  have 
slightly  lower  second  stage  ignition  times  relative 
to  CRV  shocks  -  currently  attributed  to  the  infiu- 
ence  of  the  first  stage  ignition  pressure  increase. 

Finally,  compare  the  experimental  data  to  the 
simulations.  At  high  temperatures,  the  CRV  and 
conventional  data  are  described  well  by  both  the 
UV  and  HP  computations.  This  reflects  the  accu¬ 
racy  of  both  the  LENT  mechanism  [7]  and  the 
suitability  of  simulating  single-stage  ignition  delay 
times  using  the  UV  constraint  when  the  times  are 
relatively  short  and  non-ideal  shock  tube  effects 
and  pre-ignition  energy  release  are  minimal. 
(Nonetheless,  we  regard  HP  modeling  as  physi¬ 
cally  more  correct  in  reactive  experiments  where 
the  measured  pressure  is  constant.)  At  tempera¬ 
tures  below  725  K,  it  appears  that  the  CRV  and 
conventional  complete  ignition  data  are  more  clo¬ 
sely  modeled  by  the  HP  simulation  than  by  the 
UV  modeling.  In  contrast,  above  725  K,  there  is 
slight  evidence  that  the  UV  simulation  may  more 
closely  match  the  conventional  data,  while  the  HP 
model  more  closely  tracks  the  CRV  points.  The 
difference  between  the  measured  second  stage  igni¬ 
tion  delay  times  and  the  UV  simulations,  at  tem¬ 
peratures  below  725  K,  can  once  again  be  linked 
to  the  pressure  profiles  of  the  two  experiment 
types.  At  these  low  temperatures,  the  ELNE  UV 
simulation  predicts  a  first  stage  ignition  pressure 
increase  of  as  much  as  40%,  whereas  that  observed 
in  the  CRV  experiments  was  negligible  and  that  in 
the  conventional  experiments  was  typically  only 
10%  at  all  temperatures,  refiecting  the  reality  that 
refiected  shock  ignition  is  really  not  a  constant 
UV  process.  In  the  regime  above  725  K,  the  UV 
model  predicts  a  first  stage  ignition  pressure 
increase  of  only  about  10%,  which  more  closely 
matches  that  observed  in  the  conventional  experi¬ 
ments,  whereas  the  CRV  data  were  taken  in  con¬ 
stant  pressure  conditions.  Finally,  observe  that 
the  EENE  mechanism  predicts  a  slightly  higher 
temperature  (775  K)  than  the  measurements 
(765  K)  at  which  the  minimum  first  stage  ignition 
time  occurs.  Further  experiments,  including  spe¬ 
cies  and  temperature  time-histories,  are  needed 
to  fully  characterize  this  interesting  NTC  region. 

3.2.  Fuel  time-histories  and  first  stage  fuel  decom¬ 
position  fraction 

Additional  information  on  the  details  of  igni¬ 
tion  can  be  obtained  by  measuring  time-histories 
of  important  reactants,  radicals,  intermediates, 
and  products  during  oxidation.  In  the  present 
work,  we  report  our  first  measurements  of  fuel, 
OH,  aldehydes,  CO2,  and  H2O,  as  well  as  temper¬ 
ature  time-histories,  all  conducted  in  the  CRV 
mode  of  operation.  These  measurements  should 
provide  stronger  constraints  on  the  kinetic  mech¬ 
anism,  allowing  improvements  to  be  made  with 
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Fig.  3.  HeNe  absorbance  profiles,  normalized  to  unity 
at  time  zero,  for  several  temperatures  (CRV  operation). 
Conditions:  «-heptane/15%02/5%C02/Ar,  (/)  =  0.74, 

6.5  atm.  Open  circles  denote  first  stage  ignition  (defined 
as  the  point  of  maximum  slope),  solid  triangles  denote 
post-ignition  fuel  absorbance,  and  solid  circles  denote 
second  stage  ignition. 


more  precision  than  by  simply  matching  global 
validation  targets,  e.g.  ignition  delay  times. 

HeNe  laser  absorbance  time-histories  are  plot¬ 
ted  at  several  different  NTC  region  temperatures 
in  Fig.  3.  Profiles  have  been  normalized  to  unity 
at  time  zero  in  order  to  facilitate  comparison. 
Prior  to  first  stage  ignition,  some  of  these  mea¬ 
surements  have  high  uncertainty  due  to  high 
absorbance  values  (around  a  =  4  in  some  cases). 
However,  after  first  stage  ignition,  absorbance 
values  have  fallen  due  to  fuel  decomposition  such 
that  accurate  fuel  mole  fraction  values  can  be 
deduced  after  accounting  for  interfering  species 
(which  also  absorb  at  3.39  pm)  and  for  the  effect 
of  first  stage  ignition  temperature  rise  on  the 
absorption  cross  section 

Several  observations  can  be  made  from  Fig.  3. 
First,  the  traces  show  that  first  stage  ignition  delay 
time  decreases  with  increasing  T  until  about 
765  K,  at  which  point  first  stage  ignition  delay 
time  increases.  A  similar  pattern  can  be  deduced 
for  second  stage  ignition,  although  some  shocks 
showed  scatter  in  ignition  time,  as  discussed 
above.  Furthermore,  though  convoluted  by  tem¬ 
perature  change  and  interference,  it  appears  that 
the  fraction  of  fuel  remaining  after  first  stage  igni¬ 
tion  increases  with  increasing  temperature. 

In  order  to  clarify  this  latter  observation,  an 
analysis  was  performed  to  quantify  the  fraction 
of  fuel  remaining  after  first  stage  ignition.  Using 
the  LLNL  mechanism  [7]  with  HP  modeling, 
intermediate  hydrocarbon  mole  fractions,  as  well 
as  the  temperature  and  volume  increase  upon  first 
stage  ignition,  were  estimated.  Also,  using  cross 
section  values  measured  by  Klingbeil  et  al.  [16], 
high-temperature  absorbance  values  of  interfering 
species  were  specified  and  used  to  correct  the 


observed  HeNe  absorbance.  Values  for  the 
fraction  of  fuel  remaining  after  first  stage  ignition 
have  been  plotted  in  Fig.  4.  As  also  noticed  in  [5], 
the  fraction  of  fuel  remaining  after  first  stage 
ignition  increases  with  temperature.  Both  the  data 
and  the  model  follow  an  exponential  trend,  as 
shown  in  the  semilog  scale  of  this  plot.  At  high 
temperatures,  the  LLNL  mechanism  agrees  with 
the  values  measured  in  this  work,  while  it 
under-predicts  these  values  at  low  NTC  region 
temperatures. 

A  sensitivity  analysis  performed  using  the 
LLNL  mechanism  [7]  with  HP  constraints  at 
700  K  and  7  atm  (0  =  0.75)  revealed  that  dur¬ 
ing  first  stage  ignition  the  /i-heptane  mole  frac¬ 
tion  is  subject  to  strong  competition  between 
ROO'^QOOH  rearrangement  reactions  (promot¬ 
ing  ignition)  and  ROO>^Q  +  HO2  reactions 
(inhibiting  ignition).  Further  analysis  is  necessary 
to  determine  the  extent  to  which  this  competition 
determines  the  fraction  of  fuel  remaining  after 
first  stage  ignition.  The  analysis  also  revealed  that 
the  fuel  mole  fraction  during  first  stage  ignition  is 
sensitive  to  the  n-heptane  +  OH>^R  +  H2O  and 
QOOHOO>^OQOO  +  OH  reactions. 

3.3.  OH  and  aldehydelketenelformic  acid  time- 
histories 

Absorbance  profiles  of  OH  and  aldehyde/ 
ketene/formic  acid  were  measured  during  /i-hep- 
tane  oxidation  at  initial  temperatures  of  669  and 
705  K  in  n-heptane/15%02/5%C02/Ar  mixtures 
with  an  equivalence  ratio  of  0.72  using  the  CRV 
strategy.  To  quantify  these  measurements,  the 
LLNL  mechanism  was  used  to  predict  the  post¬ 
ignition  temperatures,  enabling  specification  of 
the  proper  absorption  cross  section  values.  The 


Fig.  4.  Fraction  of  fuel  remaining  after  first  stage 
ignition  (points),  along  with  comparison  with  HP 
simulation  using  LLNL  mechanism  [7]  (dashed  line) 
and  best-fit  line  to  data  (solid  line).  Conditions: 
«-heptane/15%02/5%C02/Ar,  </>  =  0.75,  6.5  atm. 
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results  at  705  K  are  shown  in  Fig.  5,  together  with 
HP  simulations  by  the  LLNL  mechanism  [7].  In 
this  figure,  first  stage  ignition  is  marked  by  a  rapid 
rise  in  aldehyde  concentration,  and  second  stage 
ignition  is  marked  by  the  rapid  rise  in  OH  concen¬ 
tration  (to  a  predicted  plateau  level  for  OH  that  is 
dictated  by  equilibrium  thermodynamics)  and  a 
characteristically  rapid  consumption  of  intermedi¬ 
ates.  Notice  that  at  these  conditions  the  LLNL 
mechanism  under-predicts  these  times,  consistent 
with  the  under-prediction  of  ignition  times  shown 
in  Fig.  2.  However,  the  measured  plateau  values 
of  aldehyde  and  OH  mole  fraction  match  those 
of  the  model  to  within  10%. 

3.4.  CO 2  and  H2O  time-histories 

Mole  fractions  of  carbon  dioxide  and  water 
were  measured  during  ^-heptane  oxidation  at  ini¬ 
tial  temperatures  of  699  and  744  K  in  n-heptane/ 
15%02/Ar  mixtures  (no  seeded  CO2)  with  an 
equivalence  ratio  of  0.74  using  the  CRV  strategy. 
The  results  at  699  K  are  shown  in  Fig.  6,  together 
with  HP  simulations  by  the  LLNL  mechanism  [7]. 
First,  notice  that  the  first  stage  ignition  delay  time, 
as  indicated  by  the  rapidly  rising  H2O  and  CO2 
mole  fractions  around  9  ms,  is  under-predicted 
by  the  mechanism,  consistent  with  our  other  mea¬ 
surements.  Moreover,  the  complete  ignition  time, 
indicated  by  the  second  CO2  step  increase,  is  also 
under-predicted  (Unfortunately,  strong  IR  emis¬ 
sion  obscured  the  second-stage  ignition  measure¬ 
ment  of  water.)  Note  that  the  CO2  plateau 
values  following  first  and  second  stage  ignition 
closely  match  those  predicted  by  the  mechanism, 
while  the  measured  initial  H2O  plateau  value  is 
about  16%  lower  than  the  LLNL  prediction. 


Fig.  5.  Measurements  of  OH  and  aldehyde/ketene/ 
formic  acid  mole  fraction  at  an  initial  (pre-first  stage 
ignition)  temperature  of  705  K  (solid  lines),  together 
with  HP  simulation  by  the  LLNL  mechanism  [7]  (dotted 
lines).  Conditions:  «-heptane/15%02/5%C02/Ar,  (p  = 
0.72,  6.6  atm. 


Fig.  6.  Mole  fraction  of  carbon  dioxide  and  water  at  an 
initial  (pre-first  stage  ignition)  temperature  of  699  K 
(solid  lines),  together  with  HP  simulation  by  the  LLNL 
mechanism  [7]  (dotted  lines).  Conditions:  /z-heptane/ 
15%02/Ar,  (P  =  0.74,  7.0  atm. 


Sensitivity  analysis  performed  for  this  experiment 
revealed  that  the  mole  fraction  of  CO2  is  sensitive 
to  a  competition  between  ROO>^QOOH  and 
n-heptane  +  OH<^R  +  H2O  reactions,  but  further 
analysis  is  needed  to  quantitatively  determine  the 
effect  of  this  competition  on  the  post-first  stage 
ignition  plateau  level  of  CO2. 

3.5.  Temperature  time-histories  and  first  stage 
temperature  change 

First  stage  temperature  measurements  were 
carried  out  using  the  water  produced  in  the  exper¬ 
iments.  (Unfortunately,  strong  IR  emission  in 
these  experiments  obscured  the  second-stage  igni¬ 
tion  event.)  These  temperature  traces  are  shown  in 
Fig.  7,  together  with  HP-constrained  simulation 
with  the  LLNL  mechanism.  Several  points  are 
evident.  First,  observe  that  the  temperature  mea¬ 
sured  after  first  stage  ignition  very  closely  matches 
that  of  the  simulation.  Second,  despite  beginning 
at  a  lower  initial  temperature,  the  post-first  stage 
ignition  temperature  of  the  699  K  experiment  is 
essentially  identical  to  that  of  the  744  K  experi¬ 
ment;  this  is  shown  in  both  the  simulation  and 
the  experiment.  Finally,  notice  that  the  measured 
temperature  remains  relatively  constant  after  first 
stage  ignition,  whereas  that  predicted  by  the  HP 
modeling  increases  steadily  with  time,  presumably 
due  to  incorrect  parameters  of  the  mechanism. 

Using  the  thermodynamically  calculated 
pre-gnition  temperature  as  a  baseline,  the  temper¬ 
ature  increase  during  first  stage  ignition  was  com¬ 
puted  for  the  two  temperature  measurements 
discussed  above.  The  results,  together  with  HP 
predictions  by  the  LLNL  mechanism  [7],  are 
shown  in  the  inset  of  Fig.  7.  Interestingly,  the  sim¬ 
ulation  predicts  an  essentially  linear  relationship 
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Fig.  7.  Temperature  time-histories  at  initial  conditions 
(pre-first  stage  ignition)  of  699  K/7.0  atm  (dashes)  and 
744  K/6.9  atm  (solid),  together  with  HP  simulation 
using  the  LLNL  mechanism  [7].  Mixture:  w-heptane/ 
15%02/Ar,  =  0.74.  Inset:  temperature  increase  during 

first  stage  ignition  as  a  function  of  initial  (pre-first  stage 
ignition)  temperature  (points),  together  with  HP  simu¬ 
lation  with  the  LLNL  mechanism  [7]  (dashed  line). 

in  the  temperature  increase  with  the  initial  pre- 
first  stage  ignition  temperature.  The  data  match 
very  closely  to  the  simulation,  giving  confidence 
in  the  LLNL  mechanism-based  temperature-com¬ 
pensation  technique  used  to  determine  absorption 
cross  section  values  after  the  first  ignition  event  in 
the  species  time-history  analyses  above.  Although 
only  limited  temperature  data  were  observed  in 
this  study,  the  value  of  such  data  is  seen  to  be 
high,  and  hence  we  anticipate  refining  this  method 
and  employing  it  more  widely  in  future  work. 


indistinguishable,  and  the  LLNL  mechanism 
consistently  under-predicts  the  measured  times. 
Significant  differences  were  seen  between  HP  and 
UV  modeling  of  second  stage  ignition  times,  with 
HP  simulations  generally  in  better  agreement  with 
data.  The  multi-species  and  temperature  capabili¬ 
ties  demonstrated  in  this  study,  together  with  the 
staged-driver  filling  technique  enabling  long  test 
times  and  the  CRV  strategy  allowing  improved 
gasdynamic  modeling,  should  prove  very  useful 
in  future  studies  of  combustion  chemistry. 
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Appendix  A.  Supplementary  data 

Supplementary  data  associated  with  this  article 
can  be  found,  in  the  online  version,  at  http:// 
dx.doi.org/10.1016/j.proci.2014.05.001. 


4.  Summary 

We  report  here  the  first  constant-pressure 
^-heptane  ignition  delay  time  data  taken  at  6.5- 
7.4  atm  in  the  NTC  region  using  a  second-genera¬ 
tion  CRV  technique  combined  with  staged-driver 
gas  filling.  Using  these  techniques,  we  have 
achieved  test  times  of  up  to  55  ms  and  have 
observed  both  first  and  second  stage  ignition  phe¬ 
nomena  under  near-constant-pressure  conditions. 
Comparisons  between  the  conventional  and 
CRV  filling  strategies  reveal  that  ignition  delay 
data  taken  using  these  strategies,  for  the  current 
reactive  mixtures  and  test  conditions,  agree  rea¬ 
sonably  well  within  the  scatter  of  the  measure¬ 
ments.  In  addition  to  ignition  delay  times,  we 
measured  time-histories  of  fuel,  OH,  aldehydes, 
CO2,  H2O,  and  temperature  in  the  NTC  region. 
At  high  temperatures,  we  found  that  the  LLNL 
mechanism  [7]  with  either  UV  or  HP  modeling 
agrees  well  with  the  measured  single  stage  ignition 
delay  data.  In  the  NTC  region,  two  stages  of 
ignition  were  observed.  Constant  UV  and  HP 
simulations  for  first  stage  ignition  times  are 
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Ignition  delay  times  were  measured  behind  reflected  shock  waves  in  a  shock  tube  for  a  wide  variety  of 
distillate  fuels  over  a  range  of  temperatures,  pressures  and  mixtures.  The  fuels  studied  include:  jet  fuels 
(JP-5,  JP-8,  and  Jet  A),  rocket  propellants  (RP-2),  diesel  fuels  (F-76  and  DF-2)  and  gasoline.  A  simple  cor¬ 
relation  was  found  to  describe  the  ignition  delay  times  for  all  these  fuel/air  experiments  for  equivalence 
ratios  near  unity,  temperatures  from  1000  to  1400  K,  and  pressures  from  6  to  60  atm.  A  simple  correlation 
was  also  found  for  low-fuel-concentration  experiments  diluted  in  argon.  Previously  published  ignition 
delay  time  data  were  found  to  be  in  good  agreement  with  these  correlations.  Finally,  for  several  fuels 
studied,  systematic  variations  were  seen  in  the  activation  energy  for  ignition  delay  time  measurements 
with  varying  equivalence  ratio  and  oxygen  concentration. 

©  2016  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

The  quality  and  composition  of  distillate  fuels  throughout  the 
world  varies  significantly.  With  this  variation  comes  a  concern  that 
engine  performance  will  be  critically  affected  by  fuels  with  widely 
varying  composition.  For  certain  engines  such  as  gas  turbines,  sig¬ 
nificant  changes  in  ignition  delay  times,  for  example,  can  have  a 
significant  impact  on  engine  operation.  Similarly,  high- 
performance  spark  ignition  engines  are  sensitive  to  fuel  composi¬ 
tion  and  ignition  delay  time.  For  other  engines,  such  as  rockets, 
the  assumption  that  “mixed  is  burnt”  means  that  for  the  bulk  com¬ 
bustion  volume,  chemical  kinetics  are  not  assumed  to  play  a  criti¬ 
cal  role.  Chemical  kinetics  can  however  play  an  important  role  in 
predicting  rocket  engine  oscillations  and  stability.  In  military  engi¬ 
nes,  the  “One  Fuel”  directive  implemented  by  the  U.S.  Department 
of  Defense,  requires  an  understanding  of  the  influence  of  using  this 
one  fuel,  i.e.  jet  fuel,  in  engines  designed  for  other  fuels,  such  as 
diesel,  further  motivating  the  need  to  understand  variations  in 
ignition  delay  time  across  a  range  of  distillate  fuels. 

Thus,  there  continues  to  be  a  wide  research  effort  to  develop 
accurate  models  for  the  combustion  kinetics  of  jet  fuels  and  other 
distillate  fuels  with  particular  interest  on  the  role  and  importance 
of  fuel  composition.  Critical  to  this  development  is  a  need  for  accu¬ 
rate  kinetics  target  data  to  validate  and  refine  these  models.  Com¬ 
mon  laboratory  modeling  targets  include  ignition  delay  time  (IDT) 
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and  flame  speed,  which  place  global  constraints  on  model  predic¬ 
tions,  and  species  time-histories  that  place  more  specific  con¬ 
straints  on  the  internal  sub-mechanisms  of  these  models. 
Although  a  wide  literature  on  ignition  delay  time  measurements 
in  shock  tubes  exists,  these  data  are  primarily  for  neat  hydrocar¬ 
bons  and  to  a  much  lesser  extent  for  a  range  of  practical  distillate 
fuels.  For  a  review  of  methods  and  data,  for  example,  see  Lifshitz 
[1]  and  Davidson  and  Hanson  [2]. 

One  concern  of  model  developers  is  that  because  distillate  fuels 
have  variable  compositions,  kinetics  models  and  fuel  surrogates 
may  have  to  be  tuned  to  the  specific  composition  of  the  individual 
fuel  to  accurately  reproduce  critical  combustion  parameter  such  as 
ignition  delay  times  and  flame  speeds.  Based  on  our  current  knowl¬ 
edge,  for  example  of  the  large  variation  of  ignition  delay  times  for 
archetypal  individual  surrogate  components,  this  concern  is  justi¬ 
fied  [3]. 

Distillate  fuels,  however,  are  not  composed  of  a  small  and  lim¬ 
ited  number  of  components  as  are  surrogate  fuels.  Surrogate  fuels 
are  often  designed  to  include  a  small  number  of  selected  individual 
archetypal  components  of  n-alkane,  iso-alkane,  cyclo-alkane  and 
aromatic  species.  GCxGC  analyses  of  distillate  fuels,  on  the  other 
hand,  demonstrate  an  immense  array  of  different  component  spe¬ 
cies.  While  this  large  variability  in  composition  might  be  viewed 
with  concern,  recent  measurements  of  distillate  fuel  pyrolysis 
products  have  demonstrated  that  only  a  small  number  of  common 
intermediate  species  form  [4],  providing  the  basis  for  a  possible 
argument  that  some  combustion  characteristics  of  distillate  fuels 
should  actually  be  similar. 
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Related  to  this,  recent  work  on  model  development  for  practical 
distillate  fuels  (e.g.  the  hybrid  model  of  Wang  et  al.  [5])  assumes 
that  ignition  may  be  viewed  at  high  temperatures,  simply  yet  accu¬ 
rately,  as  initial  pyrolysis  (including  abstraction  reactions)  fol¬ 
lowed  by  oxidation  of  stable  intermediates.  This  assumption  is 
completely  consistent  with  shock  tube  measurements  of  the  spe¬ 
cies  time-history  of  fuel  and  intermediate  fuel  product  species 
(e.g.  the  major  product  C2H4)  during  n-heptane,  n-dodecane  and 
n-hexadecane  oxidation  [6-8].  With  this  model,  the  observed  com¬ 
monality  in  decomposition  products  may  work  in  the  modeler’s 
favor  to  simplify  the  task  of  generating  a  common  mechanism 
for  all  jet  fuels.  In  short,  if  most  fuels  generate  approximately  the 
same  decomposition  products,  should  not  the  ignition  delay  times 
be  similar? 

This  question  has  been  difficult  to  study  because  of  the  paucity 
of  high-temperature  IDT  data  for  distillate  fuels.  In  part,  this  is 
because  of  a  lack  of  a  single  source  for  research  fuel  and  the  inher¬ 
ent  variability  of  the  generic  fuel  composition.  Gauthier  et  al.  mea¬ 
sured  ignition  delay  times  in  gasoline  in  2004  [9].  Early  jet  fuel 
studies  were  performed  by  Dean  et  al.  [10].  Vasu  et  al.  in  our  lab¬ 
oratory  then  investigated  IDT  for  a  series  of  jet  fuels  [11];  this  work 
was  continued  by  Zhu  et  al.  for  F-76  [12],  and  Haylett  et  al.  for  DF-2 
[13].  Wang  and  Oehlschlaeger  [14]  and  Dooley  et  al.  [15]  provided 
high-pressure  ignition  delay  time  for  a  specific  Jet  A  fuel 
(POSF4658).  Liang  et  al.  measured  ignition  delay  times  for  China 
#3  aviation  kerosene  [16];  Steil  et  al.  measured  highly  dilute  kero¬ 
sene  in  02/argon  [17].  Other  workers  have  studied  jet  fuel  ignition 
delay  times,  but  without  a  major  emphasis  on  gas-phase  high  tem¬ 
perature  IDT  (see  Ref.  [11]). 

Recently,  Gowdagiri  et  al.  found  in  a  study  of  F-76  and  bio-fuels 
that  high-molecular  weight,  mostly  aliphatic  fuels  had  very  similar 
high-temperature  ignition  properties  [18].  In  their  study,  they 
found  an  IDT  correlation  for  these  fuels  at  temperatures  above 
1 000  K  of  the  form 

Tign  =  1.7  X  exp(25.82/R[kcal/mol  iqi)  [s]  (1) 

This  expression  captures  the  weak  dependence  on  equivalence 
ratio  ((|))  and  an  approximate  expected  dependence  on  pressure  P 
in  atm  and  temperature  T  in  K  for  these  distillate  fuels. 

In  this  study,  we  further  investigate  the  observation  of  Gowda¬ 
giri  et  al.  about  high-molecular-weight  fuels  and  measure  the  igni¬ 
tion  delay  times  for  a  wide  variety  of  distillate  fuels  over  a  wide 
range  of  temperatures,  pressures,  and  mixtures.  We  find  that  under 
certain  constraints,  these  ignition  delay  times  can  all  be  well- 
correlated  with  simple  relationships. 

2.  Experimental  method 

Ignition  delay  times  were  measured  behind  reflected  shock 
waves  in  two  shock  tubes.  High-pressure  measurements  (greater 
than  10  atm)  were  made  in  a  heated  5.0  cm  inner  diameter  shock 
tube  using  scribed  aluminum  diaphragms;  low-pressure  measure¬ 
ments  (near  1.2  atm)  were  made  in  a  heated  14.1  cm  inner  diame¬ 
ter  shock  tube  using  polycarbonate  diaphragms.  The  experimental 
setup  for  the  high-pressure  shock  tube  is  shown  in  Fig.  1.  Both 
shock  tubes  are  helium  driven,  turbo-molecular  pumped,  and 
heated  to  90  °C;  their  accompanying  gas  mixing  assemblies  were 
also  heated  to  90  °C  and  turbo-molecular  pumped.  Full  evaporation 
of  the  fuel  components  was  verified  visually  by  first  evaporating 
small  liquid  samples  of  the  fuel  in  an  external  glass  volume  at 
the  same  temperature  and  then  expanding  this  gas  into  the  mixing 
assembly  at  an  even  lower  pressure.  Separate  tests  were  performed 
to  ensure  that  the  partial  pressure  in  the  mixing  assembly  scaled 
linearly  with  the  liquid  sample  volume  over  the  range  of  fuel  sam¬ 
ple  volumes  used  for  each  experiment  (typically  0.5-2  cc).  Mano- 


metric  measurements  (in  the  mixing  assembly)  and  laser 
absorption  measurements  (in  the  shock  tube)  of  the  fuel  loading 
were  then  compared  to  ensure  that  condensation  losses  in  the 
transfer  and  shock  tube  filling  processes  were  negligible.  Incident 
shock  speeds  in  both  shock  tubes  were  measured  using  arrival 
times  from  a  series  of  5  PCB™  pressure  transducers  approximately 
equally-spaced  and  located  near  the  end  wall.  Reflected-shock 
temperatures  and  pressures  based  on  the  measured  shock  speed 
extrapolated  to  the  end  wall  and  the  ideal  frozen-chemistry  shock 
relations,  assuming  vibrationally-equilibrated  test  gas  in  both  the 
incident  and  reflected  regions,  were  accurately  determined 
(2a  =  1%).  The  assumption  that  the  reflected  shock  conditions  are 
vibrationally  equilibrated  was  tested  and  confirmed  by  Gauthier 
et  al.  in  heptane/air  mixtures  [9].  Non-reactive  pressure  profiles 
in  these  experiments  (using  N2  instead  of  O2  in  the  mixture)  were 
adjusted  using  driver  inserts  to  limit  pressure  (and  hence  also  tem¬ 
perature)  variations  to  less  than  1%  over  the  required  test  times 
[19].  (A  representative  non-reactive  pressure  profile  is  shown  in 
Fig.  2.) 

Ignition  delay  times  were  measured  using  OH*  emission  near 
306  nm  (Thorlabs™  model  PDA36A  photodiode,  UG5  Schott™  glass 
Alter,  and  simple  collection  optics)  and  side-wall  pressure  using  a 
Kistler™  603B  PZT.  Ignition  delay  times  were  defined  as  the  time 
from  the  arrival  of  the  reflected  shock  wave  at  the  side-wall  obser¬ 
vation  port  (2  cm  from  the  end  wall  for  the  low-pressure  shock 
tube  and  1  cm  from  the  end  wall  for  the  high-pressure  shock  tube) 
to  the  time  derived  from  back  extrapolating  the  steepest  rise  in  the 
OH*  or  side-wall  pressure  signal  to  the  local  baseline.  Ignition 
delay  times  using  the  two  methods  were  normally  in  close  agree¬ 
ment  (within  ±3%)  with  each  other. 

The  largest  uncertainty  in  the  IDT  measurements  can  be  associ¬ 
ated  with  the  influence  of  the  reflected  shock  temperature  on  the 
ignition  delay  times.  With  activation  energies  for  the  majority  of 
the  experiments  performed  being  near  40  kcal/mol,  uncertainties 
of  1%  in  T5,  the  reflected  shock  temperature,  translates  into  an 
absolute  uncertainty  in  the  ignition  delay  times  of  ±15%.  (Error 
bars  are  not  shown  on  the  figures  as  they  are  approximately  the 
same  dimensions  as  the  symbols.)  However,  systematic  differences 
in  the  data  of  10%  (e.g.  because  of  pressure  or  fuel/oxygen  loading) 
can  readily  be  distinguished.  Representative  data  for  JP-8  and  Jet  A 
are  shown  in  Fig.  2.  In  the  JP-8  example,  the  ramp  in  3.39  pm 
absorbance  seen  in  the  incident  shock  regime  from  -50  to  0  ps  is 
a  result  of  a  small  transverse  misalignment  of  this  beam  across 
the  shock  tube.  The  temporally-broad  pressure  rise  seen  from 
-10  to  30  ps  is  a  result  of  a  bifurcation  of  the  reflected  shock  wave 
traveling  through  the  boundary  layer  that  is  seen  with  diatomic 
gases.  The  pressure  rise  seen  near  200  ms  before  the  ignition  event 
at  210  ms  is  a  result  of  weak  non-local  ignition  increasing  the  local 
pressure.  Individual  experiments  were  examined  for  non-local 
ignition  occurring  away  from  the  end  wall  and  deflagration-to- 
detonation  transitions  (DDT)  that  can  shorten  measured  IDT.  Mea¬ 
surements  where  there  is  evidence  of  significant  shortening  of  the 
measured  IDT  were  not  included  in  the  dataset.  IDT  can  also  be 
estimated  from  rates  of  fuel  consumption  measured  using  IR  laser 
absorption  (see  next  paragraph).  Measurements  where  the  effects 
of  non-local  ignition  or  DDT  start  to  become  evident  at  the  same 
time  (within  ~5%)  as  the  estimated  IDT  (from  IR  laser  absorption) 
are  deemed  acceptable. 

In  order  to  minimize  uncertainty  in  fuel  concentration,  this 
quantity  was  measured  directly  in  the  shock  tube  using  IR  laser 
absorption  at  3.39  pm  and  Beer’s  law;  absorption  coefficients  for 
each  fuel  that  were  measured  in  separate  experiments  in  our  lab¬ 
oratory  using  the  method  of  Klingbeil  et  al.  [20].  Typical  uncertain¬ 
ties  in  the  fuel  loading  based  on  these  absorption  coefficients  have 
a  Standard  Error  of  2a  =  2-5%.  For  all  fuels  tested,  CnHm  values, 
needed  for  accurately  determining  gas-phase  equivalence  ratios 
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Fig.  2.  Representative  ignition  delay  time  data.  Upper  panel:  JP-8  (Al)/air,  initial 
reflected  shock  conditions:  1208  K,  10.2  atm,  cj)  =  1.0;  lower  panel:  Jet  A  (A2)/4%  O2/ 
argon,  1278  K,  15.3  atm,  cj)  =  1.0. 


1400  K,  and  for  cj)  near  unity,  little  sensitivity  to  (|)  was  seen  in 
the  IDT  data.  See  Section  3.7  for  a  further  discussion  of  this  point. 


3.  Results  and  discussion 

The  correlations  of  the  data  are  based  on  recent  ignition  delay 
time  measurements  in  a  variety  of  fuels.  A  summary  of  these  fuels 
and  their  properties  is  given  in  Table  1.  Table  2  outlines  the  com¬ 
positions  of  the  first  5  fuels.  A  table  of  measured  ignition  delay 
times  can  be  found  in  the  supplementary  material  associated  with 
this  paper. 

When  analyzing  the  data  and  evaluating  expressions  for  corre¬ 
lations,  the  data  was  binned  into  several  groups,  i.e.  of  high  pres¬ 
sure  (HP),  low  pressure  (LP),  Air  versus  02/Argon,  variable 
equivalence  ratio  (cj)),  and  variable  O2  mole  fraction  (X02). 

New  HP  air  IDT  data:  distillate  fuels  in  air  at  high  pressures  (10- 
60  atm)  with  d)  =  1  (0.85-1.15). 

Previous  HP  air  IDT  data:  previously  published  IDT  data  for  dis¬ 
tillate  fuels  in  air  at  high  pressures  (10-60  atm)  with  d)  =  l 
(0.85-1.15). 

HP  4%  02lAr  IDT  data:  distillate  fuels  in  4%  O2  in  argon  at  high 
pressure  (10-60  atm)  with  0  =  1  (0.85-1.15). 

LP  4%  02/Ar  IDT  data:  distillate  fuels  in  4%  O2  in  argon  at  low 
pressure  (1-15  atm)  with  0  =  1  (0.85-1.15). 

HP  N2  vs  Ar  IDT  data:  comparison  of  IDT  of  high  Pressure 
(12  atm)  distillate  fuels  in  air  and  21%02/argon. 

Variable  X02  IDT  data:  comparison  of  IDT  with  varying  X02,  O2 
concentration. 

Variable  0  IDT  data  comparison  of  IDT  with  varying  cf),  equiva¬ 
lence  ratio. 

3.1.  New  HP  air  IDT  data 

All  the  new  ignition  delay  time  (IDT)  data  in  synthetic  air  for 
distillate  fuels  (A1  through  GIO  in  Table  1)  for  equivalence  ratios 
between  0.85  and  1.15  and  for  temperatures  above  1000  K  and 
pressures  from  10  to  60  atm  can  be  correlated  with  the  simple 
formula 

Tign  =  1.037  X  exp(30.54/RT)  (2) 


in  the  mixtures,  were  determined  from  independent  H/C  ratio  and 
GCxGC  composition  measurements  [21-23]. 

Finally,  while  subtle  changes  in  IDT  are  expected  with  changes 
in  equivalence  ratio  (cf)),  experiments  were  binned  to  include  all 
measurements  with  equivalence  ratios  between  0.85  and  1.15.  This 
is  because  in  the  temperature  range  of  the  present  study,  1000- 


where  Xign  is  in  ps,  P  is  in  atm  and  the  activation  energy  is  in  kcal/- 
mol  and  temperature  are  in  Kelvins.  Correlations  were  determined 
by  a  linear  regression  of  the  logarithm  of  the  IDT  (all  equally 
weighted)  vs  the  logarithms  of  the  pressure  and  inverse  tempera¬ 
ture.  Ignition  delay  times  the  same  units  are  also  used  in  Eqs.  (3) 
and  (4).  Experiments  were  binned  to  include  equivalence  ratios 
between  0.85  and  1.15,  as  over  this  temperature  range  and  equiva¬ 
lence  ratio,  little  sensitivity  to  x  was  seen  in  the  IDT  data  (see  Sec- 
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Table  1 

Fuels  for  which  IDT  were  measured  or  cited  in  this  study. 


FUEL/identifier 

P  range  [atm] 

Equivalent  molecular  formula 

Military  POSF#  (when  available) 

Lower  heating  value  (MJ/kg) 

Aromatic%  (v/v) 

Reference 

New  data 

JP-8  ( Al ) 

1.4-35.9 

C10.8H21.8 

10,264 

43.1 

13.4 

Present  study 

Jet  A  (A2) 

0.9-14.5 

C11.4H22.1 

10,325 

43.0 

18.7 

Present  study 

JP-5  (A3) 

1.0-32.3 

C11.9H22.6 

10,289 

43.0 

20.6 

Present  study 

RP-2  (R4) 

1.0-52.2 

C12.0H24.1 

7688 

43.6 

0.3 

Present  study 

RP-2  (R5) 

1.1-15.9 

C12.6H25.7 

5433 

43.8 

0.6 

Present  study 

Kerosene  (K6) 

12.5-59.3 

C11.5H25.0 

Present  study 

Kerosene  (K7) 

21.4-61.6 

C12.0H23.0 

Present  study 

DF-2  (D8) 

54.6-59.3 

C13.7H23.0 

Present  study 

DF-2  (D9) 

11.7-56.9 

C13.6H25.9 

12,407 

42.9 

19.3 

Present  study 

Gasoline  (GIO) 

13.1-27.8 

C7.4H13.7 

42.9 

30.7 

Present  study 

Older  data 

JP-8  (Jll) 

5.6-37.8 

C10.9H22.0 

6169 

16 

[11] 

JetA(J12) 

22.0-50.5 

C10.17H19.91 

4658 

[11] 

JP-8  (J13) 

18.4-30.7 

C11H21  est. 

13.9 

[11] 

JetA(J14) 

21.8-26.2 

C11H22  est. 

[11] 

F-76  (FI 5) 

15.9-44.0 

C14.8H26.9 

[11] 

DF-2  (D16) 

4.0-7.9 

C12H23  est. 

[13] 

Gasoline  (G17) 

15.4-54.6 

C7.5H13.9  est. 

30  est. 

[9] 

JetA(J18) 

8.1-9.5 

C11H22  est. 

21.5 

[10] 

JetA(J19) 

16.7-24.8 

C10.17H19.91 

4658 

[14,15] 

Kerosene  (K20) 

5.0-27.6 

C9.49H19.5 

7.0 

[16] 

F-76  (F21) 

10.0-21.3 

C14.0H25.1 

16 

[18] 

Table  2 

Composition  breakdown  for  selected  fuels.  Composition  shown  in  volume  %.  Values  adapted  from  Refs.  [23,24]. 

Aromatics 

iso-paraffins 

n-paraffins 

cyclo-paraffins 

Al  10264 

13.4 

39.7 

26.8 

20.0 

A2  10325 

18.7 

29.5 

20.5 

31.9 

A3  10289 

20.6 

18.1 

13.9 

47.4 

R4  7688 

0.3 

38.3 

2.6 

58.3 

R5  5433 

0.6 

41.4 

14.2 

43.7 

tion  3.7  for  a  discussion  of  this  point).  The  Standard  Error  for  this 
correlation  a  =  18%  and  for  the  activation  energy  is  480  K.  The 
new  data  correlated  with  this  expression  are  shown  in  Fig.  3.  Of  par¬ 
ticular  note  is  that  the  IDT  for  the  kerosene  fuels  (A1-K8),  gasoline 
(GIO)  and  diesel  fuel  (D9)  exemplars  all  correlate  well  with  this 
expression.  These  fuels  span  molecular  formulae  from  C7.4H13.7 
(GIO)  to  C13.6H25.9  (D9)  and  contain  between  0.3%  (R4)  and  30.7% 
(GIO)  aromatic  components.  This  insensitivity  to  aromatic  compo- 


1250K  1111K  1000K 


Fig.  3.  Ignition  delay  time  data  for  distillate  fuel/air,  cj)  =  0.85-1.15,  normalized  to 
12  atm  using  the  pressure  dependence  of  Eq.  (2).  Eq.  (2)  correlation  is  shown  as  a 
solid  line. 


nents  is  consistent  with  measurements  of  heptane/toluene  mix¬ 
tures  performed  by  Schulz  et  al.  who  found  little  change  in  IDT 
from  the  neat  heptane  values  until  a  significant  fraction  of  the  mix¬ 
ture  was  aromatic  [24].  Uncertainties  in  the  mole  fraction  of  heavy 
aromatic  fuel  components,  as  a  result  of  their  high  sensitivity  to 
condensation  losses  in  the  mixing,  transfer,  and  shock  tube  filling 
processes,  are  thus  not  expected  to  have  a  significant  influence  on 
IDT  measurements. 

3.2.  Previous  HP  air  IDT  data 

Previous  distillate  fuel/air  shock  tube  ignition  delay  time  data 
(J11-J14,  J18,  J19,  G17)  from  our  laboratory  are  compared  with  this 
correlation  (Eq.  (2))  in  Fig.  4.  These  data  exhibit  greater  scatter,  but 
are  still  in  reasonable  agreement  with  Eq.  (2).  Notably,  earlier  mea¬ 
surements  of  the  same  fuel  (J19)  by  Wang  and  Oehlschlaeger 
[14,15]  and  (J12)  by  Vasu  et  al.  [9]  are  in  good  agreement  with  this 
correlation.  Good  agreement  is  also  seen  with  the  F-76  data  (F21) 
from  Gowdagiri  et  al.  [18].  The  largest  deviation  is  seen  in  the 
shortest-time  ignition  delay  measurements  (J18)  by  Dean  et  al. 
[10]  that  occurred  at  the  highest  temperatures  (above  1200  K);  lar¬ 
ger  uncertainties  are  expected  at  the  shorter  IDT  times  reported  in 
that  study. 

3.3.  HP  4%  02lAr  IDT  data 

Fig.  5  includes  high  pressure  (10-16  atm  and  50-55  atm)  distil¬ 
late  fuel  data  measured  in  4%  02/balance  argon  (Al,  A2,  R4,  R5,  D9, 
Jll,  El  5).  These  data  (with  equivalence  ratios  between  0.85  and 
1.15)  can  be  correlated  with  the  similar  formula 
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Fig.  4.  Ignition  delay  time  data  for  previously  published  fuel/air  data,  cj)  =  0.85-1 .1 5, 
normalized  to  12  atm  using  pressure  dependence  of  Eq.  (2).  Eq.  (2)  correlation  is 
shown  as  a  solid  line. 


1428K  1250K  1111K  1000K 


Fig.  5.  Ignition  delay  time  data  for  fuel/4%  02/argon  data,  cj)  =  0.85-1.15,  normalized 
to  12  atm  using  pressure  scaling  of  Eq.  (3).  Eq.  (3)  correlation  is  shown  as  a  solid 
line.  Earlier  data  are  shown  with  open  symbols. 

Tign  =  1.387  X  10-^P^’  ^°=exp(34.81/RT)  (3) 

Notably,  the  ignition  delay  time  data  in  argon  exhibit  smaller  scat¬ 
ter  (Standard  Error  is  14%)  as  shock  wave  experiments  in  a  mona¬ 
tomic  carrier  gas  do  not  experience  large  bifurcation  in  the 
reflected  shock  regime  or  suffer  significantly  from  possible  vibra¬ 
tional  relaxation  effects.  Also  notable  is  that  the  activation  energy 
(Standard  Error  =  0.68  kcal/mol)  of  the  lower  concentration  4%  O2/ 
balance  argon  experiments  is  slightly  larger  than  the  activation 
energy  for  the  fuel/air  data.  Both  sets  of  earlier  data  (Jll  &  El 5) 
are  in  excellent  agreement  with  the  correlation  [11,12]. 

3.4.  LP  4%  02lAr  IDT  data 

Eig.  6  includes  low  pressure,  1-10  atm,  distillate  fuel  data  (AI¬ 
RS,  D9)  measured  in  4%  O2  balance  argon.  These  data,  with  temper¬ 
atures  between  1000  and  1400  K  and  equivalence  ratios  between 
0.85  and  1.15,  can  be  correlated  with  the  formula 

Tign  =  4.73  X  io-3p-°=®=  exp(33.40/RT)  (4) 


Fig.  6.  Ignition  delay  time  data  for  fuel/4%  02/argon  data,  f  =  0.85-1.15,  normalized 
to  12  atm  using  pressure  dependence  of  Eq.  (4).  Eq.  (4)  correlation  is  shown  as  a 
solid  line.  (Data  has  been  scaled  to  12  atm  for  graphic  comparison  with  other  data 
in  this  study.) 

Again,  the  activation  energy  (Standard  error  is  0.81  kcal/mol)  is 
slightly  higher  for  this  lower  concentration  (lower  fuel  and  oxygen 
and  lower  pressure)  data.  At  temperatures  above  1400  K,  a  higher 
activation  energy  (52.07  kcal/mol)  would  better  capture  behavior 
of  data  for  the  jet  fuels  shown;  see  the  dashed  line  in  Eig.  6.  It  is 
important  to  note  that  while  the  scaled  IDT  values  in  this  segment 
of  the  plot  are  less  than  100  ps,  the  actual  measured  values  at  lower 
pressures  are  longer  than  100  ps,  and  the  uncertainties  in  the  IDT 
values  are  not  different  than  the  lower  temperature  data.  An 
increase  in  the  reactivity  and  shorter  IDT  values  at  the  higher  tem¬ 
peratures  is  likely  related  to  the  decreased  competition  to  radical 
formation  and  branching  (that  is  directly  linked  to  the  rate  of  the 
reaction  of  H  +  O2  =  OH  +  0)  from  radical  removal  or  scavenging 
processes  (that  are  dominated  by  reactions  of  pyrolysis  and  oxida¬ 
tion  intermediates  with  H  and  OH).  These  radical  removal  processes 
are  strongly  sensitive  to  the  distribution  of  fuel  decomposition 
products  and  to  their  decomposition  products,  both  of  which  are 
expected  to  change  with  elevated  temperature.  Also  of  note  is  that 
for  this  low-pressure  data,  the  pressure  scaling  is  significantly 
reduced  from  the  high-pressure  region. 

Because  of  the  weak  dependence  on  equivalence  ratio  for  rich 
mixtures  (discussed  in  Section  3.7),  data  at  higher  equivalence 
ratio,  up  to  (f)  =  2,  were  also  found  to  be  in  agreement  with  this 
correlation. 

3.5.  HP  N2  vs  Ar  IDT  data 

Eig.  7  shows  a  comparison  of  ignition  delay  times  at  high  pres¬ 
sures  (near  15  atm)  for  Al,  A2  and  A3  fuels  with  21%  O2  in  nitrogen 
(synthetic  air)  or  21%  O2  in  argon  as  the  carrier  gas,  along  with  the 
Eq.  (2)  correlation.  Because  of  the  lower  heat  capacity  (Cp)  for  the 
monatomic  carrier  gas  compared  to  the  diatomic  nitrogen  carrier 
gas,  temperatures  modeled  using  either  a  constant  P  or  constant 
V  assumption  can  change  differently  and  in  different  directions 
depending  on  whether  the  reaction  chemistry  is  endothermic  or 
exothermic.  In  the  case  of  fuel  ignition,  during  fuel  decomposition, 
an  endothermic  process,  the  temperature  drop  with  argon  as  the 
carrier  gas  is  larger  than  in  nitrogen.  However,  during  the  oxida¬ 
tion  process  which  is  exothermic,  the  temperature  rise  is  larger 
in  argon  than  in  nitrogen.  These  two  effects  have  an  opposite  effect 
on  the  overall  ignition  delay  time.  Argon  lengthens  the  decomposi¬ 
tion  step,  but  shortens  the  oxidation  step.  The  same  changes  hap- 
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Fig.  7.  A  comparison  of  ignition  delay  time  data  for  fuel/2 1%  02/argon  and  fuel/air 
normalized  to  12  atm. 

pen  with  nitrogen,  but  to  a  lesser  degree.  However,  independent  of 
these  possible  effects,  the  measured  values  of  the  ignition  delay 
times  in  21%02/argon  are  very  close  to  the  IDT  values  in  air,  though 
there  is  some  evidence  that  the  A2  and  A3  fuel/2 1%  02/argon  are 
tending  to  be  shorter  than  the  air  data  at  lower  temperatures. 


3.6.  Variable  X02  IDT  data 

Fig.  8  shows  a  comparison  of  IDT  values  for  jet  fuel  (A2)  with 
variable  O2  concentration  for  equivalence  ratios  near  unity.  The 
variation  in  O2  concentration  scaling  is  temperature  dependent. 
The  activation  energy  appears  to  drop  as  the  fuel  and  oxygen  con¬ 
centration  are  increased.  In  the  case  of  the  1.2  atm  4%  02/argon 
experiments,  these  data  fall  in  the  quadratic  roll-off  region  seen 
in  Fig.  5. 

At  1 2  atm  in  argon,  the  ratio  of  the  ignition  delay  times  for  A2 
fuel  with  4%  and  21%  O2  (a  ratio  of  ~l/5.25)  implies  an  IDT  scaling 
of  X52'^^  in  the  overlapping  temperature  window.  A  simple  corre- 
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Fig.  8.  A  comparison  of  A2  jet  fuel  ignition  delay  time  data  for  cj)  =  1  for  varying 
amounts  of  O2  in  argon  and  N2.  The  500  ppm  fuel/02/Ar  data  is  POSF6169  (Jet  A) 
which  is  very  similar  to  the  A2  fuel. 


lation  of  IDT  with  O2  concentration  over  the  full  concentration 
set,  however,  is  not  yet  possible. 

3.7.  Variable  0  IDT  data 

The  ignition  delay  time  correlations  in  the  Sections  3.1 -3.6 
were  limited  to  experiments  with  equivalence  ratios  between 
0.85  and  1.15.  This  is  because  an  examination  of  a  larger  set  of  data 
with  a  wider  range  of  equivalence  ratios  shows  only  limited  varia¬ 
tion  of  the  IDT  with  equivalence  ratio  in  this  range.  Fig.  9  shows  the 
variation  of  ignition  delay  time  with  equivalence  ratio  for  A1  jet 
fuel  for  two  pressures:  1.2  and  12  atm. 

What  is  evident  in  Fig.  8  is  that  the  dependence  of  ignition  delay 
time  on  equivalence  ratio  varies  with  temperature.  In  addition,  this 
variation  produces  a  cross-over  point  or  temperature,  near  1200  K, 
where  there  is  effectively  no  dependence  of  ignition  delay  time  on 
equivalence  ratio  for  A2  jet  fuel.  Above  1200  K,  as  equivalence  ratio 
is  decreased,  ignition  delay  times  get  shorter.  Below  1200  K,  as 
equivalence  ratio  is  decreased,  ignition  delay  times  get  longer.  In 
either  case,  in  this  temperature  range,  small  changes  in  equiva¬ 
lence  ratio  near  values  of  unity,  do  not  strongly  affect  the  ignition 
delay  times.  In  contrast  to  this  however.  Dean  et  al.  [7]  found  in 
their  study  at  8. 1-9.5  atm  that  higher  equivalence  ratios  always 
resulted  in  shorter  IDT  values,  and  no  cross-over  temperature 
was  observed  in  their  data  up  to  1 660  K.  Further  work  is  needed 
to  characterize  this  cross-over. 


4.  Conclusions 

Ignition  delay  times  for  a  wide  variety  of  distillate  fuels  were 
found  to  be  well  correlated  if  care  was  taken  to  identify  the  impor¬ 
tant  variables  and  areas  of  commonality.  Of  importance  is  the 
observation  that  ignition  delay  times  for  distillate  petroleum  fuels 
in  air  with  stoichiometries  near  unity  scale  together  when  pressure 
and  temperature  are  considered  for  temperatures  above  1 000  K. 
Large  variations  in  fuel  composition  do  not  appear  to  significantly 
affect  the  ignition  delay  times.  This  observation  is  useful  in  sup¬ 
porting  the  approach  of  using  a  common  detailed  or  reduced 
mechanism  to  describe  high-temperature  jet  fuel  oxidation. 

Equally  interesting  is  the  observation  that  at  temperatures 
above  1000  K,  there  is  only  a  weak  dependence  of  ignition  delay 
time  on  equivalence  ratio  in  the  range  from  0.85  to  1.15.  There 
does  appear  to  be  a  cross-over  temperature  near  1200  K,  where 
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Fig.  9.  Variation  of  ignition  delay  time  with  equivalence  ratio  for  Al  Jet  fuel. 
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below  this  temperature  the  ignition  shortens  with  higher  equiva¬ 
lence  ratio  and  above  this  temperature  the  ignition  time  lengthens. 

The  activation  energy  of  the  ignition  delay  time  was  also  found 
to  be  dependent  on  oxygen  concentration  for  equivalence  ratios 
near  unity,  with  lower  concentrations  experiencing  higher  activa¬ 
tion  energies. 

Finally,  it  should  be  noted  that  the  commonality  seen  in  the  IDT 
of  distillate  fuels  still  needs  to  be  tested  over  a  wider  equivalence 
ratio  range.  As  well,  though  this  commonality  in  IDT  is  seen  at 
higher  temperatures  as  seen  in  this  study  above  lOOOK,  there  is 
growing  evidence  that  at  lower  temperatures  (and  in  the  NTC  neg¬ 
ative  temperature  coefficient  regime)  that  significant  differences  in 
IDT  are  seen  between  different  distillate  fuels. 
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Abstract  A  high-speed  OH*  chemiluminescence  imag¬ 
ing  diagnostic  was  developed  to  image  the  structure  and 
homogeneity  of  combustion  events  behind  reflected  shock 
waves  in  the  Stanford  Constrained  Reaction  Volume  Shock 
Tube.  An  intensified  high-repetition-rate  imaging  system 
was  used  to  acquire  images  of  OH*  chemiluminescence 
(near  308  nm)  through  a  fused  quartz  shock  tube  end-wall 
window  at  10-33  kHz  during  the  combustion  of  ^-heptane 
(21  %  02/Ar,  0  =  0.5).  In  general,  the  imaging  technique 
enabled  observation  of  the  main  ignition  event  in  the  core  of 
the  shock  tube  that  corresponded  to  typical  markers  of  igni¬ 
tion  (e.g.,  pressure  rise),  as  well  as  localized  ignition  near 
the  wall  that  preceded  the  main  core  ignition  event  for  some 
conditions.  Case  studies  were  performed  to  illustrate  the 
utility  of  this  novel  imaging  diagnostic.  First,  by  comparing 
localized  wall  ignition  events  to  the  core  ignition  event,  the 
temperature  homogeneity  of  the  post-reflected  shock  gas 
near  the  end- wall  was  estimated  to  be  within  0.5  %  for  the 
test  condition  presented  (J  =  1159K,P  =  0.25  MPa).  Sec¬ 
ond,  the  effect  of  a  recession  in  the  shock  tube  wall,  cre¬ 
ated  by  an  observation  window,  on  the  combustion  event 
was  visualized.  Localized  ignition  was  observed  near  the 
window,  but  this  disturbance  did  not  propagate  to  the  core 
of  the  shock  tube  before  the  main  ignition  event.  Third,  the 
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effect  of  shock  tube  cleanliness  was  investigated  by  con¬ 
ducting  tests  in  which  the  shock  tube  was  not  cleaned  for 
multiple  consecutive  runs.  For  tests  after  no  cleaning  was 
performed,  ignition  events  were  concentrated  in  the  lower 
half  of  the  shock  tube.  In  contrast,  when  the  shock  tube  was 
cleaned,  the  ignition  event  was  distributed  around  the  entire 
circumference  of  the  shock  tube;  validating  the  cleaning 
procedure. 

1  Introduction 

Shock  tubes  have  been  used  for  over  60  years  in  a  wide 
range  of  experimental  capacities,  including  the  inves¬ 
tigation  of  chemical  kinetics  over  large  spans  of  tem¬ 
perature  and  pressure  [1].  Traditionally,  insight  into  the 
temperature,  pressure  and  species  concentrations  in  the 
shock  tube  is  gained  through  diagnostics  that  are  posi¬ 
tioned  at  ports  machined  into  the  side-wall  of  the  test 
section.  For  instance,  pressure  transducers  and  windows 
for  observing  path-integrated  light  emission  or  absorp¬ 
tion  can  be  positioned  at  these  port  locations  [2].  While 
these  diagnostics  have  proven  to  be  valuable,  they  can¬ 
not  provide  spatial  (radial)  information  regarding  the 
structure  and  homogeneity  of  the  test  gas.  A  particular 
need  not  met  by  traditional  diagnostics  is  to  understand 
the  facility-specific  uniformity  of  the  reflected  shock 
reaction  volume  for  exothermic  (combusting)  mixtures 
to  allow  valid  comparison  of  experimental  results  with 
chemical  kinetics  models  that  assume  homogeneous 
conditions.  There  are  many  potential  causes  of  non-uni¬ 
formity  within  the  shock  tube,  including  non-ideal  fluid 
flow,  reflected  shock  bifurcation  and  wall  heat  transfer, 
though  these  effects  are  generally  neglected  in  modeling 
reflected  shock  ignition. 
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Fig.  1  The  dimensions  of  the 
Stanford  Constrained  Reaction 
Volume  Shock  Tube  (CRVST) 
and  the  configuration  of  the 
imaging  system  are  shown 
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Previous  combustion  imaging  studies  in  shock  tubes 
have  utilized  multiple  techniques.  In  early  work,  Vermeer 
et  al.  [3]  used  schlieren  to  image  the  ignition  process  of 
hydrocarbons  near  the  end-wall  in  a  rectangular  shock 
tube.  Fieweger  et  al.  [4]  used  shadowgraph  to  visualize  the 
auto-ignition  of  several  fuels  in  a  square  shock  tube  with 
transparent  side- walls  next  to  the  end- wall.  Herzler  et  al. 

[5]  collected  a  single  visible-emission  image  shortly  after 
the  ignition  event  of  propane  through  a  small  window  in 
the  center  of  the  end-wall  of  a  shock  tube.  Penyazkov  et  al. 

[6]  mounted  two  (D  =  8  mm)  quartz  rods  into  the  end- wall 
to  collect  emission  (OH,  CH  and  C2  radicals)  data  simul¬ 
taneously  in  the  boundary  layer  and  the  center  of  a  shock 
tube  for  ethylene  auto-ignition.  More  recently,  Heufer  and 
Olivier  [7]  performed  schlieren  imaging  in  a  rectangular 
shock  tube  to  investigate  biofuel  combustion. 

Other  previous  imaging  techniques  of  note  investigated 
various  phenomena  in  shock  tubes  and  similar  experimental 
facilities.  Yoo  et  al.  [8]  performed  toluene  PLIF  through  the 
side- wall  of  a  shock  tube  using  a  transparent  and  square  test 
section  attachment  to  visualize  the  flow  fleld  and  measure 
the  temperature  fleld  in  non-reactive  mixtures  at  relatively 
low  temperatures.  Wegener  et  al.  [9]  obtained  holographic 
recordings  of  the  diaphragm  rupturing  in  an  expansion  tube 
through  optical  ports  in  the  wall.  In  related  work,  Mansfleld 
et  al.  [10]  used  a  transparent  end- wall  in  a  rapid  compres¬ 
sion  facility  to  image  the  ignition  process  and  determine 
the  existence  of  both  inhomogeneous  and  homogeneous 
auto-ignition.  These  imaging  studies  are  akin  to  imaging 
in  an  optically  accessible  engine,  where  an  understanding 
of  the  engine  operation  can  be  garnered  from  watching 
the  combustion  process  unfold  inside  the  engine  cylinder. 
For  example,  Stojkovic  et  al.  [11]  implemented  high-speed 
OH*  chemiluminescence  in  an  engine  through  a  window 
located  in  the  piston  to  study  soot  formation  and  oxidation. 

In  the  current  work,  OH*  chemiluminescence,  a  com¬ 
mon  emission  diagnostic  for  combustion  systems,  was  used 
as  a  qualitative  marker  of  the  reaction  zone.  OH*  chemi¬ 
luminescence  (peak  emission  occurring  near  308  nm)  was 
used  to  indicate  the  locations  of  high  radical  concentrations 
[12,  13].  With  emission  in  the  UV,  OH*  chemilumines¬ 
cence  was  easily  studied  with  very  little  ambient  light  con¬ 
tamination  by  using  a  band-pass  Alter. 


Here,  combustion  events  in  a  shock  tube  were  imaged 
through  a  UV-transparent  (fused  quartz)  end-wall  with  a 
high-repetition-rate  intensifled  camera.  The  resulting  col¬ 
lected  images  were  axially  integrated  emission  indicating 
the  spatial  location  and  geometry  of  the  ignition  event.  In 
this  work,  the  reaction  volume  was  constrained  to  a  small 
portion  of  the  end- wall  (see  Sect.  2.1  for  details),  which 
drastically  limited  the  length  of  the  integration  path.  For 
potentially  long  axial  integration  paths  (e.g.,  conventional 
shock  tubes),  greater  care  would  need  to  be  taken.  In  gen¬ 
eral,  chemiluminescence  images  must  be  interpreted  with 
caution  because  they  are  integrated  measurements,  but  they 
yield  powerful  information  otherwise  unavailable  from  tra¬ 
ditional  diagnostics  through  side- wall  ports  [14]. 

Through  this  work,  an  improved  understanding  of  shock 
tube  performance  and  diagnostics  was  sought.  Three  case 
studies  that  illustrate  the  utility  of  this  imaging  diagnostic 
were  performed:  Temperature  Homogeneity  (Sect.  3.1), 
Shock  Tube  Wall  Defects  (Sect.  3.2),  and  Shock  Tube 
Cleanliness  (Sect.  3.3). 

2  End- wall  imaging  system  design 

2.1  Stanford  Constrained  Reaction  Volume  Shock  T\ibe 

The  Stanford  Constrained  Reaction  Volume  Shock  Tube 
(CRVST),  as  shown  schematically  in  Fig.  1,  was  devel¬ 
oped  to  study  the  kinetics  of  low-vapor-pressure  fuels 
[15].  The  internal  diameter  of  the  buffer  and  test  sections 
is  11.53  cm,  and  these  sections  contribute  to  9.7  m  of  the 
shock  tube  length.  The  driver  section,  with  a  length  of 
3.63  m,  has  an  internal  diameter  of  17.9  cm.  A  driver  exten¬ 
sion  piece  increases  the  driver  section  length  to  a  maximum 
of  13.42  m  and  has  an  internal  diameter  of  15.4  cm  (further 
detailed  information  regarding  the  CRVST,  the  driver  sec¬ 
tion  and  driver  inserts  can  be  found  in  [16]).  Located  2  cm 
away  from  the  end-wall,  machined  ports  are  situated  that 
house  windows  for  laser  or  emission  diagnostics  and  other 
sensors  such  as  pressure  transducers. 

The  shock  tube  was  run  in  constrained  reaction  volume 
(CRV)  mode.  CRV  mode  employed  a  sliding  gate  valve  to 
mechanically  conflne  the  reactive  gas  mixture  to  a  small 
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Fig.  2  The  four-pieee  transpar¬ 
ent  end-wall  window  design 
allows  spatial  imaging  of  the 
ignition  event  behind  a  reflected 
shock  wave  through  the  end- 
wall  of  the  Stanford  CRVST 


fraction  of  the  shock  tube  length,  near  the  end- wall  [15]. 
This  gate  valve  was  opened  just  before  the  shock  wave  was 
initiated,  i.e.,  when  the  Lexan  diaphragm  burst,  in  order 
to  minimize  the  dilution  of  the  test  gas  mixture.  The  axial 
extent  of  the  test  gas  after  the  reflected  shock,  L5,  ranged 
from  5  to  10  cm  for  the  conditions  studied  [17]. 

For  each  experiment,  the  driver  section  was  filled  with 
a  tailored  helium  and  nitrogen  mixture.  The  buffer  section 
contained  a  tailored  mixture  of  carbon  dioxide,  nitrogen 
and  argon  (see  Refs.  [15,  16]  for  driver  and  buffer  gas  tai¬ 
loring  details).  For  all  tests  presented  here,  the  test  section 
was  filled  with  normal  heptane,  oxygen  and  argon,  at  an 
equivalence  ratio  of  0  =  0.5.  Research-grade  helium,  nitro¬ 
gen,  carbon  dioxide,  argon  and  oxygen  were  obtained  from 
Praxair,  Inc.,  and  spectroscopic-grade  normal  heptane  from 
Sigma-Aldrich  was  used.  Test  gas  mixture  compositions 
were  determined  via  partial  pressures  in  a  mixing  tank. 
Best  practices  previously  developed  to  minimize  the  effect 
of  mixing  of  the  test  gas  and  buffer  gas  were  followed  [17]. 

Initial  test  gas  conditions  (behind  the  reflected  shock 
wave)  were  calculated  from  the  speed  of  the  incident  shock 
wave,  determined  using  the  time-of-arrival  technique  via 
seven  side-wall  pressure  transducers  (PCB-113A26)  dis¬ 
tributed  over  the  last  section  of  the  shock  tube.  Time  zero 
for  the  experiment  was  defined  at  the  half-rise  of  a  Kistler 
(603B1)  pressure  transducer  signal  for  the  pressure  jump 
due  to  the  reflected  shock  wave.  This  pressure  transducer 
was  located  2  cm  from  the  end-wall  and  was  covered  in 
room-temperature  vulcanizing  silicone  (RTV).  For  refer¬ 
ence,  the  time  difference  between  the  incident  shock  arrival 
at  the  end- wall  and  the  passage  of  the  reflected  shock  by 
the  pressure  transducer  at  2  cm  was  typically  70|xs  and 


varies  slightly  with  initial  conditions.  Additionally,  a  tradi¬ 
tional  OH*  emission  diagnostic  at  2  cm  from  the  end-wall 
was  implemented  using  a  focusing  lens,  slit,  UG-5  Schott 
Glass  band-pass  filter  and  a  silicon  photodetector  (Thor- 
Labs  PDA36A).  Signals  from  the  pressure  transducers  and 
emission  detector  were  acquired  at  10  MHz  using  a  Lab- 
VIEW  data  acquisition  system. 

Initial  reflected  shock  temperatures  for  this  work  ranged 
from  1057  to  1240  K.  Initial  reflected  shock  pressures 
ranged  from  0.25  to  0.81  MPa. 

2.2  End-wall  window 

The  transparent  end- wall  window,  as  shown  in  Fig.  2,  con¬ 
sists  of  four  pieces:  shim,  window,  housing  and  end  cap. 
The  window  is  a  13.97  cm  (5.5  in.)  diameter,  2.54  cm 
(1  in.)  thick  cylinder.  In  this  work,  a  fused  quartz  window 
(GM  Associates)  transmits  the  308  nm  OH*  chemilumines¬ 
cence.  This  easily  removable  window  is  a  critical  design 
feature,  as  it  allows  cleaning  of  the  window  with  acetone 
between  experiments.  The  design  also  does  not  require 
any  RTV  or  sealant,  further  aiding  in  the  ease  of  window 
removal  and  re-installation.  The  window  is  rated  for  a  max¬ 
imum  pressure  of  2.4  MPa  (safety  factor  of  4).  An  impor¬ 
tant  feature  of  this  end-wall  window  design  is  that  it  allows 
optical  access  out  to  the  edge  of  the  shock  tube  wall. 

The  Teflon  shim  is  0.127  mm  (0.005  in.)  thick  and  pre¬ 
vents  over- tightening  of  the  shock  tube  end  cap.  It  allows 
the  shock  tube  to  be  adequately  sealed,  while  preventing 
excessive  compressive  stress  on  the  window  that  can  cause 
cracks.  The  2.54  mm  (0.1  in.)  thick  Teflon  housing  holds 
the  window  and  protects  the  back  of  the  window  from 
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contact  with  the  metal  of  the  shock  tube  end  cap.  The  end 
cap  interior  is  hollow  and  is  attached  to  the  test  section  of 
the  shock  tube  using  four  bolts. 

2.3  Imaging  system 

A  Vision  Research  Phantom  v710  high-speed  CMOS  cam¬ 
era  coupled  to  an  intensifier  was  used  to  collect  the  OH* 
chemiluminescence.  The  intensifier  was  a  LaVision  HS- 
IRO  equipped  with  a  gen  II  S20  photocathode.  A  Sodem 
UV  lens,  with  a  focal  length  of  /  =  100  mm,  was  used 
with  the  aperture  set  to  f/4  for  all  cases.  To  filter  out  stray 
light  and  other  emission,  an  Asahi  Spectra  high-trans¬ 
mission  band-pass  filter  centered  at  313  nm  with  10  nm 
FWHM  was  used. 

In  the  interest  of  protecting  the  camera  and  intensifier  in 
the  event  of  a  window  or  end  cap  failure,  a  UV-enhanced 
mirror  was  used  to  position  the  camera  at  a  90°  angle  to 
the  shock  tube  (Fig.  1).  Additionally,  a  shield  was  placed 
between  the  shock  tube  and  camera. 

For  the  v710  camera,  there  is  a  trade-off  between  image 
resolution  (height  x  width  pixels)  and  acquisition  rate;  a 
smaller  window  of  pixels  enables  a  faster  acquisition  rate. 
To  allow  the  full  shock  tube  cross  section  to  be  imaged  at 
a  faster  rate,  the  camera  was  moved  farther  away  from  the 
shock  tube.  In  this  study,  resolutions  were  used  that  range 
from  448  x  376  to  880  x  800  pixels,  corresponding  to 
acquisition  rates  of  33-10  kHz  (inter-frame  times  ranging 
from  0.03  to  0.1  ms),  respectively. 

For  the  highest  resolution,  the  projected  pixel  size 
was  6.59  pixels/nnn  (0.15  mm/pixel).  For  the  lowest 
resolution,  the  projected  pixel  size  was  3.26  pixels/mm 
(0.31mm/pixel).  The  depth  of  field  for  the  10  kHz  setup 
was  around  13  cm  and  for  the  33  kHz  setup  was  around 
110  cm,  both  captured  the  full  depth  of  the  compressed  test 
gas  section  (L5  =  5  cm).  The  system  was  focused  at  the 
location  of  the  diagnostic  ports,  2  cm  from  the  end-wall. 

The  intensifier  was  gated  to  4  or  5  |xs.  The  intensifier 
gain  was  adjusted  specifically  to  visualize  the  locations  of 
the  earliest  flames  with  25-50  %  of  the  maximum  possible 
camera  signal  level.  In  this  work,  the  gain  was  typically  set 
at  50-60  %  of  the  maximum  gain  setting. 

3  Results  and  case  studies 
3.1  Temperature  homogeneity 

In  a  typical  shock  tube  auto-ignition  experiment,  it  is 
assumed  that  the  test  gas  is  uniformly  heated  by  the 
reflected  shock  wave  and  rapid  reaction  occurs  homogene¬ 
ously  throughout  the  reactor  volume  at  the  ignition  delay 
time  (IDT).  The  collected  Kistler  pressure  transducer  (2  cm 


Fig.  3  a  The  pressure  trace  from  the  Kistler  pressure  transducer  and 
voltage  from  the  OH*  emission  detector  (2  cm  from  the  end-wall) 
of  the  combustion  of  n-heptane  (21  %  02/Ar,  0  =  0.5)  are  shown. 
The  initial  reflected  shock  conditions  were  1159  K  and  0.25  MPa. 
b  The  normalized  spatially  averaged  signal  from  the  collected  end- 
wall  images  is  shown.  The  images  were  recorded  at  12  kHz  with  5  p.s 
exposure  time.  The  solid  vertical  bars  were  numbered  corresponding 
to  the  timing  of  the  images  shown  in  Fig.  4 

1159K  0.25MPa 


t=  1.34ms  t=  1.43ms  t=  1.51  ms 
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Fig.  4  The  combustion  of  n-heptane  (21  %  02/Ar,  0  =  0.5)  with  ini¬ 
tial  reflected  shock  conditions  of  1159  K  and  0.25  MPa  was  recorded 
at  12  kHz  with  5  |jis  exposure  time.  A  white  outline  was  included  to 
mark  the  location  of  the  shock  tube  wall.  The  timing  of  the  three 
images  was  also  indicated  in  Fig.  3 

from  the  end- wall)  signal  is  shown  in  Fig.  3a.  The  first  rise 
in  pressure  (near  t  =  —0.1  ms)  is  due  to  the  incident  shock 
wave  passing  over  the  pressure  transducer  and  the  second 
rise  (/  =  0  ms)  is  due  to  the  reflected  shock  wave.  The  sub¬ 
sequent  large  “spike”  in  pressure  (near  /  =  1.5  ms)  corre¬ 
sponds  to  the  ignition  event  and  is  representative  of  how 
IDT  may  be  measured.  Small  changes  in  the  pressure  trans¬ 
ducer  signal  or  emission  diagnostic,  prior  to  the  dramatic 
increase  in  pressure,  are  observed  under  some  conditions 
and  may  be  taken  as  evidence  of  partial  inhomogeneity  in 
the  ignition  event  that  needs  to  be  better  understood. 

This  imaging  technique  revealed  the  spatial  distribution 
of  ignition  events  that  occured  in  the  shock  tube  as  they 
evolved  in  time.  By  acquiring  a  time  series  of  images,  the 
variability  and  hence  the  uncertainty  in  the  IDT  were  quanti¬ 
fied.  Uncertainty  in  measuring  the  IDT  is  an  often  discussed 
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topic  due  to  multiple  choices  of  definitions,  ignition  indica¬ 
tors,  and  a  large  variability  in  initial  conditions.  IDT  indica¬ 
tors  include  the  pressure  signal  or  various  species  from  laser 
absorption  diagnostics.  Definitions  of  IDT  include  the  peak, 
time  of  the  maximum  rate  of  change,  or  extrapolating  the 
maximum  slope  to  the  zero  signal  level  of  the  chosen  indi¬ 
cator  [18].  In  this  work,  the  OH*  chemiluminescence  col¬ 
lected  in  the  images  was  used  for  defining  the  IDT. 

In  the  case  presented  in  Fig.  3,  significant  (>.01  nor¬ 
malized  average  image  signal)  OH*  chemiluminescence 
was  observed  in  the  image  at  the  time  of  the  first  vertical 
line,  before  indication  of  ignition  was  observed  in  either 
the  pressure  transducer  signal  or  the  traditional  OH*  emis¬ 
sion  diagnostic  (both  located  2  cm  from  the  end- wall).  The 
CRVST  is  generally  employed  under  conditions  where 
pressure  increases  are  minimal,  thereby  facilitating  reac¬ 
tion  modeling.  Here,  higher  fuel  concentrations  and  reac¬ 
tion  temperatures  were  intentionally  chosen  to  emphasize 
inhomogeneous  effects  and  pressure  variations. 

Figure  4  shows  three  sequential  images  collected  at 
12  kHz  (inter-frame  time  of  .083  ms)  for  the  test  reported 
in  Fig.  3.  A  video  is  available  in  the  supplementary  mate¬ 
rial  that  contains  additional  images  collected  for  the  test 
depicted  in  Figs.  3  and  4. 

In  the  first  image  {t  =  1.34  ms)  of  Fig.  4,  initial  OH* 
chemiluminescence  was  observed  near  the  walls  of  the 
shock  tube  (indicated  by  the  white  outline).  Ignition  in  the 
core  was  evidenced  in  the  second  image  {t  =  1.43  ms), 
and  ignition  was  evident  everywhere  within  the  shock  tube 
cross  section  by  the  time  of  the  third  image  (f  =  1.51  ms). 

Using  the  timing  of  the  images  and  known  character¬ 
istics  of  the  combustion  event,  a  semiquantitative  esti¬ 
mate  of  the  homogeneity  within  the  shock  tube  was  made. 
Equation  1  expresses  the  exponential  temperature  depend¬ 
ence  of  the  IDT  in  an  Arrhenius-like  form.  The  relationship 
between  the  variation  (a  measure  of  uncertainty)  in  the  IDT 
and  the  variation  in  temperature  is  shown  in  Eq.  2  where 
A  is  the  pre-exponential  factor,  Ea  is  the  activation  energy, 
R  is  the  universal  gas  constant  and  T  is  the  temperature. 
Depending  on  the  IDT  correlation,  the  pre-exponential  fac¬ 
tor  can  be  expressed  as  a  function  of  multiple  experimental 
variables  (e.g.,  pressure,  oxygen  concentration  and  equiva¬ 
lence  ratio)  [19]. 


Tign  =  A  exp 


(1) 


Using  Eq.  2  and  the  assumption  that  the  pressure  was  uni¬ 
form  throughout  the  test  volume,  the  spatial  variation  in 
temperature  was  estimated  from  the  variation  of  the  IDT 
throughout  the  test  section.  The  key  assumption  here  was 


Window  Recession 


Fig.  5  Port  locations  are  equally  spaced  circumferentially  at  an  axial 
distance  of  2  cm  from  the  end- wall  in  the  Stanford  CRVST.  The 
0.76  mm  recession  caused  by  the  window  is  exaggerated  in  the  draw¬ 
ing.  The  plug  that  houses  the  flat  window  has  the  same  contour  as  the 
shock  tube  wall 


that  the  pockets  of  gas  within  the  reaction  volume  ignited 
at  different  times  due  to  the  slight  temperature  differences. 
Though  this  was  a  simplistic  assumption,  it  was  neverthe¬ 
less  useful  for  initial  quantification  of  the  temperature 
variation. 

The  variation  in  the  IDT  throughout  the  shock  tube  was 
estimated  as  the  time  interval  between  the  first  image  with 
significant  OH*  chemiluminescence  (f=  1.34  ms)  and 
the  core  ignition  event  with  the  maximum  OH*  chemilu¬ 
minescence  (f  =  1.51  ms).  This  corresponded  to  a  value 
of  Arign  =  0.17  ms.  Using  Eq.  2  and  Ea  =  186.6  kJ/mol 
(44.6  kcal/mol)  for  n-heptane,  a  value  of  AT  =  6.7  K  was 
obtained  [19].  This  estimate  indicated  that  the  temperature 
variation  throughout  the  test  section  was  around  0.5  %  for 
this  case. 

3.2  Shock  tube  wall  defects 

To  implement  laser  absorption  and  emission  diagnostics 
in  the  Stanford  CRVST,  windows  are  placed  in  machined 
ports  located  2  cm  from  the  end- wall.  There  are  eight  ports 
equally  spaced  around  the  shock  tube  circumference,  which 
can  be  filled  with  smoothly  contoured  plugs  or  plugs  con¬ 
taining  windows.  When  windows  are  installed,  a  0.76  mm 
recession  is  created  in  the  wall  at  the  top  and  bottom  of  the 
window  surface  (Fig.  5). 

These  recessions  may  cause  spatial  temperature  non¬ 
uniformities  due  to  interactions  with  the  incident  and 
reflected  shock  waves.  To  test  this  hypothesis,  window 
plugs  were  installed  at  different  circumferential  locations 
during  ignition  experiments.  Indeed,  reactions  localized 
to  the  windows  were  observed,  as  shown  in  Fig.  6.  Initial 
flames  formed  at  the  location  of  the  window,  in  addition  to 
other  locations  (some  seemed  to  correspond  to  locations  of 
smooth  plugs,  i.e.,  plugs  without  windows)  along  the  shock 
tube  wall.  A  magnified  color  scale  was  used  to  display  the 
first  image  of  the  sequence  to  highlight  the  location  of  the 
initial  flames. 
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Fig.  6  The  effects  of  the 
window  location,  evolving  over 
time,  were  recorded  at  33  kHz 
with  4  |xs  exposure  time  for  the 
combustion  of  /z -heptane  (21  % 
02/Ar,  0  =  0.5).  The  first  image 
of  each  sequence  has  a  magni¬ 
fied  color  scale  to  highlight  the 
location  of  the  initial  flames. 

The  initial  reflected  shock 
temperature  and  pressure  along 
with  the  window  location  (black 
rectangle)  were  labeled  for  both 
experiments.  A  white  outline 
was  included  to  mark  the  loca¬ 
tion  of  the  shock  tube  wall 
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These  findings  demonstrated  that  the  small  recession 
in  a  window  plug  can  be  sufficient  to  promote  initial  local 
reactions.  However,  ignition  in  the  core  seemed  unaffected. 
The  flames  that  formed  on  the  wall  did  not  propagate  to 
the  core  before  the  main  ignition  event.  Nonetheless,  it  is 
important  to  recognize  the  possibility  of  local  reactions 
near  windows  when  implementing  line-of-sight  absorption 
or  emission  diagnostics.  This  could  lead  to  a  bias  toward 
shorter  measured  ignition  delay  times  by  those  diagnostics 
and  could  be  an  explanation  for  the  early  rise  in  OH*  emis¬ 
sion  diagnostic  seen  in  Fig.  3a  for  the  previous  case  study. 

3.3  Shock  tube  cleanliness 

After  three  to  five  experiments  in  the  shock  tube,  particu¬ 
late  can  accumulate  near  the  end- wall  of  the  shock  tube. 
Procedures  for  cleaning  the  shock  tube  have  been  devel¬ 
oped  to  minimize  the  potential  effects  of  these  particulates. 
The  CRVST  was  cleaned  between  experiments  by  brushing 
out  and  cleaning  the  last  2  m  length  near  the  end- wall  with 
acetone  before  mechanically  pumping  the  shock  tube  to  a 
pressure  of  <0.1  Pa.  The  efficacy  of  the  cleaning  method¬ 
ology  was  investigated  by  imaging  the  combustion  event 
after  the  shock  tube  was  not  cleaned. 

Figure  7a  shows  two  different  combustion  events  that 
took  place  after  the  shock  tube  was  not  cleaned  for  multi¬ 
ple  runs  (the  shock  tube  was  always  pumped  down  between 
runs).  Bright  spots  were  seen  that  likely  result  from  the 
small  particulates  that  accumulated  at  the  bottom  of  the 
shock  tube.  Many  seemingly  independent  flames  were 
observed  to  initiate  in  the  bottom  of  the  shock  tube. 

Figure  7b  shows  preignition  images  from  tests  after  the 
shock  tube  was  cleaned  at  similar  initial  reflected  shock 
temperatures  and  pressures.  The  reaction  zones  were  dis¬ 
tributed  circumferentially  for  both  tests,  in  contrast  to  the 
tests  conducted  with  particulates  present. 


t=  3.74  ms  t=  3.74  ms 

1057  K  0.30  MPa  1060  K  0.30  MPa 


(b) 


Fig.  7  a  The  combustion  of  zz-heptane  (21  %  02/Ar,  0  =  0.5),  after 
the  shock  tube  was  not  cleaned,  was  recorded  at  10  kHz  with  5  p.s 
exposure  time.  The  left  image  was  acquired  after  three  experiments 
were  performed  without  cleaning,  and  the  right  image  was  acquired 
after  four  experiments  were  performed  without  cleaning,  b  The  com¬ 
bustion  of  zz-heptane  (21  %  02/Ar,  0  =  0.5),  after  the  last  2  m  of  the 
shock  tube  was  cleaned,  was  recorded  at  12  kHz  with  5  p.s  exposure 
time.  The  initial  reflected  shock  temperatures  and  pressures  were 
labeled  for  each  experiment.  A  white  outline  was  included  to  mark 
the  location  of  the  shock  tube  wall 


t  =  3.63  ms 
1093  K  0.26  MPa 


t  =  3.78  ms 
1095  K  0.26  MPa 


It  is  important  to  mention  that  for  the  other  case  stud¬ 
ies  presented  previously,  the  shock  tube  was  cleaned  before 
each  run.  Confirming  that  an  unclean  shock  tube  changes 
the  qualitative  structure  of  preignition  reactions  was  not 
surprising,  but  was  a  reminder  that  careful  control  of  exper¬ 
imental  variables  is  needed  to  achieve  homogeneous  or 
near-homogeneous  combustion  in  reflected  shock  experi¬ 
ments.  End- wall  imaging  has  the  potential  to  be  a  valuable 
tool  for  validating  other  aspects  of  shock  tube  performance 
and  for  facilitating  further  development  of  best  practices 
for  shock  tube  operation. 
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4  Conclusion 

High-speed  end-wall  imaging  of  OH*  chemilumines¬ 
cence  was  shown  to  be  a  powerful  tool  for  understanding 
time-resolved,  spatial  details  about  shock  tube  kinetics 
experiments.  The  end- wall  window  provides  new  access 
for  imaging  of  the  structure  and  location  of  the  earliest 
reaction  zones  in  the  shock  tube,  as  well  as  the  temporal 
development. 

Three  case  studies  that  illustrated  the  utility  of  the  tech¬ 
nique  were  presented.  An  estimate  of  the  temperature  vari¬ 
ability  throughout  the  test  section  was  found  to  be  around 
0.5  %  for  the  test  condition  presented  (^-heptane,  21  %  O2 
/Ar,  0  =  0.5).  The  window  recessions  were  found  to  cause 
initial  flames  to  form  at  the  window  location.  The  impor¬ 
tance  of  routine  shock  tube  cleaning  was  also  demonstrated. 

Acknowledgments  This  material  is  based  upon  work  supported 
by  the  US  Army  Researeh  Laboratory  and  the  US  Army  Researeh 
Offiee  under  Contraet/Grant  Number  W911NF1310206,  with  Ralph 
Anthenien  as  eontraet  monitor.  We  aeknowledge  loan  of  the  inten¬ 
sified  eamera  system  from  an  AFOSR-sponsored  researeh  program. 
V.A.T.  was  supported  by  the  Gabilan  Stanford  Graduate  Fellow¬ 
ship  and  the  National  Seienee  Foundation  Graduate  Researeh  Fel¬ 
lowship  Program  (Grant  No.  DGE-1 14747).  Any  opinion,  findings 
and  eonelusions  or  reeommendations  expressed  in  this  material 
are  those  of  the  authors(s)  and  do  not  neeessarily  refleet  the  views 
of  the  National  Seienee  Foundation.  M.F.C.  was  supported  by  a 
National  Defense  Seienee  and  Engineering  Graduate  Eellowship 
(32CER168a). 


References 

1.  W.  Bleakney,  D.K.  Weimer,  C.H.  Eleteher,  The  shoek  tube:  a 
faeility  for  investigations  in  fluid  dynamies.  Rev.  Sei.  Instrum. 
20(11),  807-815  (1949) 

2.  R.K.  Hanson,  D.E.  Davidson,  Reeent  advanees  in  laser  absorp¬ 
tion  and  shoek  tube  methods  for  studies  of  eombustion  ehemis- 
try.  Prog.  Energy  Combust.  Sei.  44,  103-114  (2014) 

3.  D.J.  Vermeer,  J.W.  Meyer,  A.K.  Oppenheim,  Auto-ignition  of 
hydroearbons  behind  refleeted  shoek  waves.  Combust.  Elame  18, 
327-336  (1972) 


4.  K.  Eieweger,  R.  Blumenthal,  G.  Adomeit,  Self-ignition  of  S.I. 
engine  model  fuels:  a  shoek  tube  investigation  at  high  pressure. 
Combust.  Elame  109(4),  599-619  (1997) 

5.  J.  Herzler,  L.  Jerig,  P.  Roth,  Shoek- tube  study  of  the  ignition  of 
propane  at  intermediate  temperatures  and  high  pressures.  Com¬ 
bust.  Sei.  Teehnol.  176(10),  1627-1637  (2004) 

6.  O.G.  Penyazkov,  K.L.  Sevrouk,  V.  Tangirala,  N.  Joshi,  High- 
pressure  ethylene  oxidation  behind  refleeted  shoek  waves.  Proe. 
Combust.  Inst.  32(2),  2421-2428  (2009) 

7.  K.A.  Heufer,  H.  Olivier,  S.P.  Medvedev,  S.V.  Khomik,  Optieal 
investigation  of  shoek  indueed  ignition  of  different  biofuels, 
in  Proeeedings  of  the  23rd  ICDERS,  University  of  California, 
Irvine,  CA,  24-29  July  201 1 

8.  J.  Yoo,  D.  Mitehell,  D.E.  Davidson,  R.K.  Hanson,  Planar  laser- 
indueed  fluoreseenee  imaging  in  shoek  tube  flows.  Exp.  Eluids 
49,  751-759  (2010) 

9.  M.  Wegener,  M.  Suteliffe,  R.  Morgan,  Optieal  study  of  a  light 
diaphragm  rupture  proeess  in  an  expansion  tube.  Shoek  Waves 
10,  167-178  (2000) 

10.  A.B.  Mansfield,  M.S.  Wooldridge,  H.  Di,  X.  He,  Low-tempera¬ 
ture  ignition  behavior  of  iso-oetane.  Euel  139,  79-86  (2015) 

11.  B.D.  Stojkovie,  T.D.  Eansler,  M.C.  Drake,  V.  Siek,  High-speed 
imaging  of  OH*  and  soot  temperature  and  eoneentration  in  a 
stratified-eharge  direet-injeetion  gasoline  engine.  Proe.  Combust. 
Inst.  30,  2657-2665  (2005) 

12.  V.  Nori,  J.  Seitzman,  Evaluation  of  ehemilumineseenee  as  a  eom¬ 
bustion  diagnostie  under  varying  operating  eonditions.  in  46th 
AIAA  Aerospaee  Seienees  Meeting  and  Exhibit,  Reno,  Nevada, 
7-10  Jan  2008 

13.  C.T.  Bowan,  D.J.  Seery,  Chemilumineseenee  in  the  high-temper- 
ature  oxidation  of  methane.  Combust.  Elame  12,  611-614  (1968) 

14.  H.N.  Najm,  PH.  Paul,  C.J.  Mueller,  PS.  Wyekoff,  On  the  ade- 
quaey  of  eertain  experimental  observables  as  measurements  of 
flame  burning  rate.  Combust.  Elame  113,  312-332  (1998) 

15.  M.E.  Campbell,  A.M.  Tulgestke,  D.E.  Davidson,  R.K.  Hanson,  A 
seeond-generation  eonstrained  reaction  volume  shock  tube.  Rev. 
Sei.  Instrum.  85,  055108  (2014) 

16.  M.E.  Campbell,  T.  Parise,  A.M.  Tulgestke,  R.M.  Spearrin,  D.E. 
Davidson,  R.K.  Hanson,  Strategies  for  obtaining  long  constant- 
pressure  test  times  in  shock  tubes.  Shock  Waves  25(6),  651-665 
(2015) 

17.  M.E.  Campbell,  Studies  of  biodiesel  surrogates  using  novel 
shock  tube  techniques.  Ph.D.  Thesis,  Stanford  University  (2014) 

18.  D.E.  Davidson,  R.K.  Hanson,  Interpreting  shock  tube  ignition 
data.  Int.  J.  Chem.  Kinet.  36,  510-523  (2004) 

19.  D.C.  Horning,  D.E.  Davidson,  R.K.  Hanson,  Study  of  the  high- 
temperature  autoignition  of  n-Alkane/02/Ar  mixtures.  J.  Pro¬ 
puls.  Power  18(2),  363-371  (2002) 


^  Springer 


